HU-Mathematik: 2013-06 HU-EP-13/15 
IPPP/13/09 DCPT/13/18 SL AC-PUB- 15409 
LAPTH-016/13 CERN-PH-TH/2013-058 



Leading singularities and off-shell conformal integrals 



James Drummond'^'^, Claude Duhr'^'^, Burkhard Eden^, Paul Heslop^, 
Jeffrey Pennington^, Vladimir A. Smirnov*^ '' 

" CERN, Geneva 23, Switzerland 

^ LAPTH, CNRS et Universite de Savoie, F-74941 Annecy-le-Vieux Cedex, France 

^ Institut fiir Theoretische Physik, ETH Zurich, CH-8093, Switzerland 

Institute for Particle Physics Phenomenology, 
University of Durham, Durham, DHl 3LE, U.K. 

^ Institut fiir Mathematik, Humboldt- Universitdt, Zum grofien Windkanal 6, 12489 Berlin 

^ Mathematics department, Durham University, Durham DHl 3LE, United Kingdom 

^ SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94309, USA 

^ Nuclear Physics Institute of Moscow State University, Moscow 119992, Russia 



Abstract 

The three-loop four-point function of stress-tensor multiplets in A/" = 4 super Yang-Mills 
theory contains two so far unknown, off-shell, conformal integrals, in addition to the known, 
ladder-type integrals. In this paper we evaluate the unknown integrals, thus obtaining the 
three-loop correlation function analytically. The integrals have the generic structure of 
rational functions multiplied by (multiple) polylogarithms. We use the idea of leading 
singularities to obtain the rational coefficients, the symbol - with an appropriate ansatz 
for its structure - as a means of characterising multiple polylogarithms, and the technique 
of asymptotic expansion of Feynman integrals to obtain the integrals in certain limits. 
The limiting behaviour uniquely fixes the symbols of the integrals, which we then lift 
to find the corresponding polylogarithmic functions. The final formulae are numerically 
confirmed. The techniques we develop can be applied more generally, and we illustrate this 
by analytically evaluating one of the integrals contributing to the same four-point function 
at four loops. This example shows a connection between the leading singularities and the 
entries of the symbol. 



In memory of Francis Dolan. 
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1 Introduction 

The work presented in this paper is motivated by recent progress in planar TV = 4 super 
Yang-Mills (SYM) theory in four dimensions, although the methods that we exploit and 
further develop should be of much wider applicability. 

A/" = 4 SYM theory has many striking properties due to its high degree of symmetry; 
for instance it is conformally invariant, even as a quantum theory pQ, and the spectrum of 
anomalous dimensions of composite operators can be found from an integrable system [2]. 
Most strikingly perhaps, it is related to IIB string theory on AdSsxS^ by the AdS/CFT 
correspondence |3]. This is a weak/strong coupling duality in which the same physical 
system is conveniently described by the field theory picture at weak coupling, while the 
string theory provides a way of capturing its strong coupling regime. The strong coupling 
limit of scattering amplitudes in the model has been elaborated in ref. [1] from a string 
perspective. The formulae take the form of vacuum expectation values of polygonal Wilson 
loops with light-like edges. 

This duality between amplitudes and Wilson loops remains true at weak coupling [5], 
extending to the finite terms in A/" = 4 SYM previously known relations between the 
infrared divergences of scattering amplitudes and the ultra-violet divergences of (light- 
like) Wilson lines in QCD [6]. Furthermore, it was recently discovered that both sides of 
this correspondence can be generated from n-point correlation functions of stress-tensor 
multiplets by taking a certain light-cone limit [7]. 

The four-point function of stress-tensor multiplets was intensely studied in the early 
days of the AdS/CFT duality, in the supergravity approximation [S] as well as at weak 
coupling. The one-loop [9] and two-loop [10] corrections are given by conformal ladder 
integrals. 

A Feynman-graph based three-loop result has never become available because of the 
formidable size and complexity of multi-leg multi-loop computations. Already the two 
parallel two-loop calculations [TU] drew heavily upon superconformal symmetry. However, 
a formulation on a maximal ('analytic') superspace [TT|[T2] makes it apparent that the loop 
corrections to the lowest x-space component are given by a product of a certain polynomial 
with linear combinations of conformal integrals, cf. ref. [I3HT6] . Then in ref. [I71[T8], using 
a hidden symmetry permuting integration variables and external variables, the problem of 
finding the three- loop integrand was reduced down to just four unfixed coefficients without 
any calculation and furhter down to only one overall coefficient after a little further analysis. 
This single overall coefficient can then easily be fixed e.g. by comparing to the MHV four- 
point three-loop amplitude [19] via the correlator /amplitude duality or by requiring the 
exponentiation of logarithms in a double OPE limit [T7] . 

Beyond the known ladder and the 'tennis court', the off-shell three- loop four-point 
correlator contains two unknown integrals termed 'Easy' and 'Hard' in ref. [17]. In this 
work we embark on an analytic evaluation of the Easy and Hard integrals postulating that 

• the integrals are sums Yli RiFii where Ri are rational functions and Fi are pure 
functions, i.e. Q-linear combinations of logarithms and multiple polylogarithms [20], 
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• the rational functions Ri are given by the so-called leading singularities (i.e. residues 
of global poles) of the integrals [2T] , 

• the symbol of each Fi can be pinned down by appropriate constraints and then 
integrated to a unique transcendental function. 

The principle of uniform transcendentality, innate to the planar A/" = 4 SYM theory, implies 
that the symbols of all the pure functions are tensors of uniform rank six. Our strategy will 
be to make an ansatz for the entries that can appear in the symbols of the pure functions 
and to write down the most general tensor of uniform rank six of this form. We then impose 
a set of constraints on this general tensor to pin down the symbols of the pure functions. 
First of all, the tensor needs to satisfy the integrability condition, a criterion for a general 
tensor to correspond to the symbol of a transcendental function. Next the symmetries 
of the integrals induce additional constraints, and finally we equate with single variable 
expansions corresponding to Euclidean coincidence limits. The latter were elaborated for 
the Easy and Hard integrals in ref. [22|[23] using the method of asymptotic expansion 
of Feynman integrals |2l]. This expansion technique reduces the original higher-point 
integrals to two-point integrals, albeit with high exponents of the denominator factors and 
complicated numerators. 

To be specific, up to three loops the off-shell four-point correlator is given by P|[TU|[T7] 





G4(l, 2, 3, 4) = Gf + ^^^- /^ R{1, 2, 3, 4) [aF^'^ + a^F^^) + a'F<^'^ + 0{a')] , (1.1) 



Here denotes the number of colors and a is the 't Hooft coupling. represents the 
tree-level contribution and -R(l,2,3,4) is a universal prefactor, in particular taking into 
account the different SU{4:) flavours which can appear (see ref. [171[T8] for details). Our 
focus here is on the loop corrections. These can be written in the compact form (exposing 
the hidden S^+i symmetry) as 




, X2, Xs, Xi) 




12-^13-^14-^23'^24-^34 



/ 





where 




(1.3) 




^^12^34^56 + permutations 



(1.4) 



l<i<i<6 "^ij 
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(^12)^(^34^45^56^67^73) + ^7 permutations 




(1.5) 



l<i<j<T^ij 
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Writing out the sum over permutations in the above expressions, these are written as 
follows 

i^^'^=^7l234, (1.6) 
F^^^ = hi2;34 + ^34;12 + /i23;14 + /il4;23 (1-7) 
+ ^13;24 + ^24;13 + 2 (-^12^34 + ^13^24 + ^U^ls) [^'1234] ^ , 

F^^^ = [Li2;34 + 5 perms ] + [Ti2;34 + 11 perms ] (1.8) 
+ [-E'i2;34 + 11 perms ] + l[xl^xl^Hi2-34 + 11 perms ] 
+ [(fi' X h) 12-34 + 5 perms ] , 
which involve the following integrals: 

5'1234 — ^ / 2 2 2~ ' y^-^ 

46) 



rT7-4 / f /Y«2 fy»2 ^2 "\ ^2 / ^2 fy»2 ^2 
J l"''15"''35"''45J"''56l"''26"''36"'' 



At three-loop order we encounter 

^12^34 f d'^X5d'^XQd'^X7 



{g X h) 



12;34 - R 



-77-6 / f fy»2 Ay«2 fy»2 Ay«2 "\ / fy* 2 ™2 fy» 2 "\ / ™2 ^2 fy»2 \ ^ 

~' l"''15"''25"''35"''45Jl"''16-^36'^46Jl-^27-^37"''47M' 

X34 /" (i'^Xs (i'^Xe (i'^Xy 

12;34 — / / 2 ^2 ^2 N^2 /'^2 ^2 N^2 /'^2 ^2 ^2 N ' 



67 



/' /Y«2 fy»2 ^2 "\ ^2 / ^2 fy»2 \ fy»2 / fy»2 ^2 

l'''15"''35"''45J"''56l"''36-^46J"''67l"''27'''37'''' 



47; 



rp ^ ^ f C^^Xsrf^XeC/^Xy xjj 

12,34 ^6 / ( rf'^ ^2 '\/^2 ^2 \/^2 ^2 ^2 \^2 ^2 ^2 ' 

'I J l"''15"''35Jl"''16"''46A"''37"''27"''47M56"''57"''67 



E 



12:34 



/TT-G / o^2 /y>2 ^2 ^^2 /" ^2 /Y>2 ^2 \^2 / /y>2 ^2 ,-yi2 ^ 
J l"''15"''25"''35M56l"''26"''36-^46M67l"''17"''27"''47J 

X34 /" (i^Xs d^X% d'^X-j X57 



12,34 /tt6 / /y»2 ^2 /y>2 ^2 \^2 / ^2 /y»2 \^2 / ^2 /yi2 ^2 ,yi2 

l"''15"''25"''35"''45M56l"''36"''46M67l"''17"''27"''37"''' 



47 



Here h, L are recognised as the one-loop, two-loop and three-loop ladder integrals, respec- 
tively, the dual graphs of the off-shell box, double-box and triple-box integrals. Off-shell, 
the 'tennis court' integral T can be expressed as the three- loop ladder integral L by using 
the conformal flip propertied of a two-loop ladder sub-integral [25] . The only new integrals 
are thus E and H (see fig. [T]). 



^Such identities rely on manifest conformal invariance and will be broken by the introduction of most 
regulators. For instance, the equivalence of T and L is not true for the dimensionally regulated on-shell 
integrals. 
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-£^12;34 



-f^l2;34 



Figure 1: The Easy and Hard integrals contributing to the correlator of stress tensor 
multiplets at three loops. 



Conformal four-point integrals are given by a factor carrying their conformal weight, 
say, (3:^3X24)" times some function of the two cross ratios 



.= ^ = .0:, . = |i^ = (l-.)(l-a:). (1.11) 

X]_3X24 "''13"''24 



Ladder integrals are explicitly known for any number of loops, see ref. [26] where they 
are very elegantly expressed as one-parameter integrals. Integration is simplified by the 
change of variables from the cross-ratios {u, v) to ( X, clS defined in the last equation. The 
unique rational prefactor, x1^X24 {x — x), is common to all cases and can be computed by 
the leading singularity method as we illustrate shortly. This is multiplied by pure polylog- 
arithm functions which fit with the classification of single-valued harmonic polylogarithms 
(SVHPLs) in ref. [27]. The associated symbols of the ladder integrals are then tensors 
composed of the four letters {x, x, 1 — x, 1 — x}. 

On the other hand, for generic conformal four-point integrals (of which the Easy and 
Hard integrals are the first examples) there are no explicit results. Fortunately, in recent 
years a formalism has been developed in the context of scattering amplitudes to find at least 
the rational prefactors (i.e. the leading singularities), which are given by the residues of the 
integrals [21] . There is one leading singularity for each global pole of the integrand and it is 
obtained by deforming the contour of integration to lie on a maximal torus surrounding the 
pole in question, i.e. by computing the residue at the global pole. As an illustratioio, let 
us apply this technique to the massive one- loop box integral (71234 defined in eq. (11. 9p . Its 
leading singularity is obtained by shifting the contour to encircle one of the global poles of 



^The massless box-integral (i.e. the same integral in the limit xf^_^_i -> 0) is discussed in ref. [55] in 
terms of twistor variables as the simplest example of a 'Schubert problem' in projective geometry. The 
ofF-shell case that we discuss here was also recently discussed by S. Caron-Huot (see [29j). 
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the integrand, where all four terms in the denominator vanish. To find this let us consider 
a change of coordinates from to Pi = xf^. The Jacobian for this change of variables is 

J = det (^^^ = det (-2x^5) , = det (4x^5 ■ xj^) = 16det (x-^- - x% - xjg) , 

(1-12) 

where the second identity follows by observing that det(M) = ^/det{MM^. Using this 
change of variables the massive box becomes 

^1234 = / J ■ (1-13) 

To find its leading singularity we simply compute the residue around all four poles at pi = 
(divided by 27ii). We obtain 



5-1234 -r^T , A1234 = \/det(x?)i •=i..4 = xf^xl^^x-x) (1.14) 



in full agreement with the analytic result 

Note that we do not consider explicitly a contour around the branch cut associated 
with the square root factor J in the denominator of (11.131) . Because there is no pole at 
infinity, the residue theorem guarantees that such a contour is equivalent to the one we 
already considered. On the other hand, in higher-loop examples, Jacobians from previous 
integrations cannot be discarded in this manner. In all the examples we consider, these 
Jacobians always collapse to become simple poles when evaluated on the zero loci of the 
other denominators and thereby contribute non-trivially to the leading singularity. 

The main results of this paper are the analytic evaluations of the Easy and Hard 
integrals. Due to Jacobian poles, the Easy integral has three distinct leading singularities, 
out of which only two are algebraically independent, though. The Hard integral has two 
distinct leading singularities, too. Armed with this information we then attempt to find 
the pure polylogarithmic functions multiplying these rational factors. Our main inputs for 
this are analytic expressions for the integrals in the limit x — )■ obtained from the results 
in [23]. Matching these asymptotic expressions with an ansatz for the symbol of the pure 
functions we obtain unique answers for the pure functions. 

The pure functions contributing to the Easy integral are given by SVHPLs, correspond- 
ing to a symbol with entries drawn from the set {x, 1 — x,x,l — x}. In this case there is a 
very straightforward method for obtaining the corresponding function from its asymptotics, 
by essentially lifting HPLs to SVHPLs as we explain in the next section. However, the 
SVHPLs are not capable of meeting all constraints for the pure functions contributing to 
the Hard integral, so that we need to enlarge the set of letters. A natural guess is to include 
X — X (cf. ref. [30]) since it also occurs in the rational factors, and indeed this turns out 
to be correct. Ultimately, one of the pure functions is found to have a four- letter symbol 
corresponding to SVHPLs, but the symbol of the other function contains the new letter: 
the corresponding function cannot be expressed through SVHPLs alone, but it belongs to 
a more general class of multiple polylogarithms. 
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Figure 2: The four-loop integral IiJ.23 defined in eq. f ll.lSp . 



Let us stress that the analytic evaluation of the Easy and Hard integrals completes 
the derivation of the three-loop four-point correlator of stress-tensor multiplets in = 4 
SYM. The multiple polylogarithms that we find can be numerically evaluated to very high 
precision, which paves the way for tests of future integrable system predictions for the 
four-point function, or for instance for further analyses of the operator product expansion. 

Finally, since our set of methods has allowed to obtain the analytic result for the Easy 
and Hard integrals in a relatively straightforward way (despite the fact that these are 
not at all simple to evaluate by conventional techniques) we wish to investigate whether 
this can be repeated to still higher orders. We examine a first relatively simple looking, 
but non-trivial, four-loop example from the list of integrals contributing to the four-point 
correlator at that order |18j : 



The computation of its unique leading singularity follows the same lines as at three loops. 
However, just as for the Hard integral, the alphabet {x, 1—x, x, 1—x} and the corresponding 
function space are too restrictive. Interestingly, this integral is related to the Easy integral 
by a differential equation of Laplace type. Solving this equation promotes the denominator 
factor 1 — M of the leading singularities of the Easy integral to a new entry in the symbol 
of the four-loop integral. Note that it is at least conceivable that the letter x — x arrives 
in the symbol of the Hard integral due to a similar mechanism, although admittedly not 
every integral obeys a simple differential equation. 
The paper is organised as follows: 

• In Section [21 we give definitions of the concepts introduced here: symbols, harmonic 
polylogarithms, SVHPLs, multiple polylogarithms and so on. 

• In Section [31 we comment on the asymptotic expansion of Feynman integrals. 




/ 




d'^x^d'^x^d^xid'^xs x\^X24^x\ 



(1.15) 



15-^18"''25-^26"''37-^38"''45'^46-^47-^48"''56-^67-^78 
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• In Sections H] and E] we derive the leading singularities, symbols and ultimately the 
pure functions corresponding to the Easy and Hard integrals. We also present nu- 
merical data indicating the correctness of our results. 

• In Section [TJ we perform a similar calculation for the four-loop integral, I^^^ . 

• Finally we draw some conclusions. We include several appendices collecting some 
formulae for the asymptotic expansions of the integrals and alternative ways how to 
derive the analytic results. 



2 Conformal four-point integrals and single- valued poly- 
logarithms 

The ladder-type integrals that contribute to the correlator are known. More precisely, if 
we write 



5'13;24 = 2^ 2 ^^^\u,v) , 
^13^24 

^13;24 = ^^^^^\U,V), (2.1) 



0C-\ oX' 



13-^24 



^13;24 =^^$('H 



^2 ^2 



then the functions ^^^\u, v) are given by the well-known result 

L\{L-iy.jQ V + [1 - u - V) ^ + u 

X (log^ + loge)'^~'(log^ + 21oge) (2.2) 



X — X \x — la; — 1 

where the conformal cross ratios are given by eq. (11. lip and where we defined the pure 
function 

/(^)(x, x) = Y^ ^^ulI^V-U ^"^'^"^"^^ (Li2w(x) - Li2i„,(x)) . (2.3) 

At this stage, the variables [x, x) are simply a convenient parametrisation which rationalises 
the two roots of the quadratic polynomial in the denominator of eq. (12. 2p . We note that 
X and X are complex conjugate to each other if we work in Euclidean space while they are 
both real in Minkowski signature. 

The particular combination of polylogarithms that appears in eq. (12. 2p is not random, 
but it has a particular mathematical meaning: in Euclidean space, where x and x are 
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complex conjugate to each other, the functions are single- valued functions of the 
complex variable x. In other words, the combination of polylogarithms that appears in the 
ladder integrals is such that they have no branch cuts in the complex x plane. In order to 
understand the reason for this, it is useful to look at the symbols of the ladder integrals. 



2.1 The symbol 

One possible way to define the symbol of a transcendental function is to consider its total 
differential. More precisely, if F is a function whose differential satisfies 

dF = J2F^dhsR^, (2.4) 



where the Ri are rational functions, then we can define the symbol of F recursively by [3T] 

S{F) = J2^iF^)®R^■ (2-5) 

i 

As an example, the symbols of the classical polylogarithms and the ordinary logarithms 
are given by 

S{Un{z)) = -{1 - z) z . . .(^ z and 5 ( ^ln"2 ) = 2 (g) . . . (g) ^ . (2.6) 



(n— 1) times 

In addition the symbol satisfies the following identities, 

...(g) (a - 6)®... = ...(ga(g... + ...(g6(g... , 

(±1)®... =0, (2.7) 
S{FG) =S{F)mS{G), 

where III denotes the shuffle product on tensors. Furthermore, all multiple zeta values 
are mapped to zero by the symbol map. Conversely, an arbitrary tensor 

Ci,...i^uJn ® • • • ® w»„ (2.8) 

ii,...,in 

whose entries are rational functions is the symbol of a function only if the following inte- 
gmhility condition is fulfilled, 

Q,„i„ dlogUi^ A d\ogUi^^^ (g . . . (g (g Wi^^j (g . . . (g Wi„ = , (2.9) 

for all consecutive pairs {ik,ik+i)- 

The symbol of a function also encodes information about the discontinuities of the 
function. More precisely, the singularities (i.e. the zeroes or infinities) of the first entries of a 
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symbol determine the branching points of the function, and the symbol of the discontinuity 
across the branch cut is obtained by dropping this first entry from the symbol. As an 
example, consider a function F{x) whose symbol has the form 

S{F{x)) = (ai - x) O . . . ® (a„ - x) , (2.10) 

where the Oj are independent of x. Then F{x) has a branching point at x = ai, and the 
symbol of the discontinuity across the branch cut is given by 

S [disCaiF(x)] = 27ri (a2 - x) O . . . (g) (a„ - x) . (2.11) 

If F is a Feynman integral, then the branch cuts of F are dictated by Cutkosky's rules. In 
particular, for Feynman integrals without internal masses the branch cuts extend between 
points where one of the Mandelstam invariants becomes zero or infinity. As a consequence, 
the first entries of the symbol of a Feynman integral must necessarily be Mandelstam 
invariants |32] . In the case of the four-point position space integrals we are considering in 
this paper, the first entries of the symbol must then be distances between two points, xfj 
for i, j = 1 . . . 4. Combined with conformal invariance, this implies that the first entries of 
the symbols of conformally invariant four-point functions can only be cross ratios. As an 
example, consider the symbol of the one-loop four-point function, 

S\f^^\x,x)] = u^^—^ + v^-. (2.12) 

'- ^ 1 — X X 

The first entry condition puts strong constraints on the transcendental functions that 
can contribute to a conformal four-point function. In order to understand this better let 
us consider a function whose symbol can be written in the form 

S{F) =u^Su + v(^S^ = (xx) ^Su + [(1 -x){l- x)] ® , (2.13) 

where Su and are tensors of lower rank. Let us assume we work in Euclidean space 
where x and x are complex conjugate to each other. It then follows from the previous 
discussion that F has potential branching points in the complex x plane at x G {0, 1, oo}. 
Let us compute for example the discontinuity of F around x = 0. Only the first term in 
eq. fl2.13p can give rise to a non-zero contribution, and x and x contribute with opposite 
signs. So we find 

S [disco (F)] = 27ri Su - 27ri Su = . (2.14) 

The argument for the discontinuities around x = 1 and x = oo is similar. We thus 
conclude that F is single-valued in the whole complex x plane. This observation puts 
strong constraints on the pure functions that might appear in the analytical result for a 
conformal four-point function. In particular, the ladder integrals ^^^^ are related to the 
single-valued analogues of the classical polylogarithms, 

Dn{x) ='nnJ2 log'l^l Lin-fc(a;) , (2.15) 

fc=0 
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where denotes the real part for n odd and the imaginary part otherwise and are 
the Bernoulli numbers. For example, we have 

= 4iD2(a;) . (2.16) 

2.2 Single- Valued Harmonic Polylogarithms (SVHPLs) 

For more general conformal four-point functions more general classes of polylogarithms may 
appear. The simplest extension of the classical polylogarithms are the so-called harmonic 
polylogarithms (HPLs), defined by the iterated integral^ [33] 

H{ai,...,an;x)= dt fa,{t) H{a2, ■ ■ . ,an;t) , aie{0,l}, (2.17) 
Jo 

with 

/o(x) = - and /i(x) = -^. (2.18) 

X 1 — X 

By definition, H{x) = 1 and in the case where all the at are zero, we use the special 
definition 

H{0^;x) = ^\og^x. (2.19) 

The number n of indices of a harmonic polylogarithm is called its weight. Note that the 
harmonic polylogarithms contain the classical polylogarithms as special cases, 

H{On-iA;x) = Un{x). (2.20) 

In ref. [35] it was shown that infinite classes of generalised ladder integrals can be ex- 
pressed in terms of single- valued combinations of HPLs. Single-valued analogues of HPLs 
were studied in detail in ref. [27], and an explicit construction valid for all weights was pre- 
sented. Here it suffices to say that for every harmonic polylogarithm of the form H{a; x) 
there is a function Cs{x) with essentially the same properties as the ordinary harmonic 
polylogarithms, but in addition it is single-valued in the whole complex x plane. We will 
refer to these functions as single-valued harmonic polylogarithms (SVHPLs). Explicitly, 
the functions Cs{x) can be expressed as 



^s{x) = ^ Cij H{di; x) H{dj; x) , (2.21) 



where the coefficients Cij are polynomials of multiple ( values such that all branch cuts 
cancel. 

There are two natural symmetry groups acting on the space of SVHPLs. The first 
symmetry group acts by complex conjugation, i.e., it exchanges x and x. The conformal 
four-point functions we are considering are real, and thus eigenfunctions under complex 



•^In the following we use the word harmonic polylogarithm in a restricted sense, and only allow for 
singularities at x G {0, 1} inside the iterated integrals. 
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conjugation, while the SVHPLs defined in ref. [27] in general are not. It is therefore 
convenient to diagonalise the action of this symmetry by defining 

L,{x) = I [Cix) - i-lpC,{x)] , 

^ (2.22) 
Ls{x) = - [C,{x) + {-lf\C,{x)] , 

where |a| denotes the weight of Cg_{x). Note that we have apparently doubled the number 
of functions, so not all the functions L^t^x) and Lci{x) can be independent. Indeed, one can 
observe that 

Ld{x) = [product of lower weight SVHPLs of the form Lg{x) ] . (2.23) 

The functions L^^x) can thus always be rewritten as linear combinations of products of 
SVHPLs of lower weights. In other words, the multiplicative span of the functions Lct{x) 
and multiple zeta values spans the whole algebra of SVHPLs. As an example, in this basis 
the ladder integrals take the very compact form 

/(^)(a:,x) = (-l)^+^2Lo_o,o,i,o_o(^) • (2-24) 

L-l L-1 

While we present most of our result in terms of the Lg{x), we occassionally find it convenient 
to employ the Ld{x) and the Cs{x) to obtain more compact expressions. 

The second symmetry group is the group 5*3 which acts via the transformations of the 
argument 

x^x, x^l-x, z -> 1/(1 - x) , (2.25) 
X— )-l/x, X— J-l — 1/x, X— 7'X/(x — 1). 

This action of 5*3 permutes the three singularities {0, 1, oo} in the integral representations 
of the harmonic polylogarithms. In addition, this action has also a physical interpretation. 
The different cross ratios one can form out of four points Xj are parametrised by the group 
Sa/ (Zi2 X Z2) ~ 5*3. The action (12.251) is the representation of this group on the cross ratios 
in the parametrisation (11.111) . 

2.3 The x^O limit of SVHPLs 

We will be using knowledge of the asymptotic expansions of integrals in the limit x — 
in order to constrain, and even determine, the integrals themselves. If the function lives 
in the space of SVHPLs there is a very direct and simple way to obtain the full function 
from its asymptotic expansion. 

This direct procedure relies on the close relation between the series expansion of SVH- 
PLs around x = and ordinary HPLs. In the case where SVHPLs are analytic at (x, x) = 
(i.e. when the corresponding word ends in a '1') then 

lim£^(x) = H^{x) . (2.26) 
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Similar results exist in the case where is not analytic at the origin. In that case 

the limit does strictly speaking not exist, but we can, nevertheless, represent the function 
in a neighbourhood of the origin as a polynomial in log u, whose coefficients are analytic 
functions. More precisely, using the shuffie algebra properties of SVHPLs, we have a unique 
decomposition 

(x) = ^ ttp^^' logP u C^' (x) , (2.27) 

p,w' 

where integer numbers and are analytic at the origin = 0. 

Conversely, if we are given a function /(x, x) that around x = admits the asymptotic 
expansion 

= ^ap,^log''Mi7^(x) + 0(x), (2.28) 

where the independent of (x, x) and w are words made out of the letters and 

1 ending in a 1, there is a unique function fsy}iPh{x,x) which is a linear combination of 
products of SVHPLs that has the same asymptotic expansion around x = as f{x,x). 
Moreover, this function is simply obtained by replacing the HPLs in eq. f l2.28p by their 
single-valued analogues, 

) = J]ap,«,log^M/:^(a;) . (2.29) 

p,w 

In other words, f{x,x) and /svhpl(2:,x) agree in the limit x — up to power-suppressed 
terms. 

It is often the case that we find simpler expressions by expanding out all products, i.e. 
by not explicitly writing the powers of logarithms of u. More precisely, replacing logw by 
log X -|- log X in eq. (12.281) and using the shuffie product for HPLs, we can write eq. (I2.28P 
in the form 

/(x,x) = ^a^i/^(x) + logxP(x,logx) + C(x), (2.30) 

w 

where P(x,logx) is a polynomial in logx whose coefficients are HPLs in x. From the 
previous discussion we know that there is a linear combination of SVHPLs that agrees 
with /(x, x) up to power-suppressed terms. In fact, this function is independent of the 
actual form of the polynomial P, and is completely determined by the first term in the 
left-hand side of eq. f lfgU]) . 

/svHPL(a:, x) = ^ Cn,{x) . (2.31) 

w 

So far we have only described how we can always construct a linear combination of 
SVHPLs that agrees with a given function in the limit x — )■ up to power-suppressed 
terms. The inverse is obviously not true, and we will encounter such a situation for the 
Hard integral. In such a case we need to enlarge the space of functions to include more 
general classes of multiple polylogarithms. Indeed, while SVHPLs have symbols whose 
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entries are all drawn from the set {x,x, 1 — x, 1 — x}, it was observed in ref. [30] that 
the symbols of three-mass three-point functions (which are related to conformal four-point 
functions upon sending a point to infinity) in dimensional regularisation involve functions 
whose symbols also contain the entry x — x. Function of this type cannot be expressed 
in terms of HPLs alone, but they require more general classes of multiple polylogarithms, 
defined recursively by G{x) = 1 and, 

G{ai,...,an;x)= / G(a2, • • • , a„; t) , G{Op;x) = ^-^, (2.32) 

where Oj G C. We will encounter such functions in later sections when constructing the 
analytic results for the Easy and Hard integrals. 

3 The short-distance limit 

In this section we sketch how the method of 'asymptotic expansion of Feynman integrals' 
can deliver asymptotic series for the x — )■ limit of the Easy and the Hard integral. These 
expansions contain enough information about the integrals to eventually fix ansatze for the 
full expressions. 

In ref. [22|[23] asymptotic expansions were derived for both the Easy and Hard integrals 
in the limits where one of the cross ratios, say u, tends to zero. The limit n — )■ 0, — ?■ 1 
can be described as a short-distance limit, X2 — )■ Xi. Let us assume that we have got rid 
of the coordinate X4 by sending it to infinity and that we are dealing with a function of 
three coordinates, xi,X2,X3, one of which, say xi, can be set to zero. The short-distance 
limit we are interested in then corresponds to X2 — 0, so that the coordinate X2 is small 
(soft) and the coordinate X3 is large (hard). This is understood in the Euclidean sense, i.e. 
X2 tends to zero precisely when each of its component tends to zero. One can formalise 
this by multiplying X2 by a parameter p and then considering the limit p — > upon which 
u ~ p^, — 1 ~ p. 

For a Euclidean limit in momentum space, one can apply the well-known formulae 
for the corresponding asymptotic expansion written in graph-theoretical language (see 
ref. [2^ for a review). One can also write down similar formulae in position space. In 
practice, it is often more efficient to apply the prescriptions of the strategy of expansion 
by regions p^[36] (see also Chapter 9 of ref. [37j for a recent review), which are equivalent 
to the graph-theoretical prescriptions in the case of Euclidean limits. The situation is 
even simpler in position space where we work with propagators l/x^^. It turns out that 
in order to reveal all the regions contributing to the asymptotic expansion of a position- 
space Feynman integral it is sufficient to consider each of the integration coordinates x, 
either soft (i.e. of order X2) or hard (i.e. of order X3). Ignoring vanishing contributions, 
which correspond to integrals without scale, one obtains a set of regions relevant to the 
given hmit. One can reveal this set of regions automatically, using the code described in 
refs. [MlEn]. 
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The most complicated contributions in the expansion correspond to regions where the 
internal coordinates are either all hard or soft. For the Easy and Hard integrals, this 
gives three-loop two-point integrals with numerators. In ref. [22], these integrals were 
evaluated by treating three numerators as extra propagators with negative exponents, so 
that the number of the indices in the given family of integrals was increased from nine 
to twelve. The integrals were then reduced to master integrals using integration-by-parts 
(IBP) identities using the C++ version of the code FIRE [10] • While this procedure is not 
optimal, it turned out to be sufficient for the computation in ref. [22]. In ref. [23], a more 
efficient way was chosen: performing a tensor decomposition and reducing the problem to 
evaluating integrals with nine indices by the well-known MINCER program [H], which is 
very fast because it is based on a hand solution of the IBP relations for this specific family 
of integrals. This strategy has given the possibility to evaluate much more terms of the 
asymptotic expansion. 

It turns out that the expansion we consider includes, within dimensional regularisation, 
the variable u raised to powers involving an amount proportional to e = (4 — (i)/2. A 
characteristic feature of asymptotic expansions is that individual contributions may exhibit 
poles. Since the conformal integrals we are dealing with are finite in four dimensions, the 
poles necessarily cancel, leaving behind some logarithms. The resulting expansions contain 
powers and logarithms of u times polynomials in f — 1. Instead of the variable v, we turn 
to the variables {x,x) defined in eq. (11. lip . Note that it is easy to see that in terms of 
these variables the limit u ^ 0,v ^ 1 corresponds to both x and x becoming small. 

In fact, we only need the leading power term with respect to u and all the terms with 
respect to x. The results of ref. [23] were presented in terms of infinite sums involving 
harmonic numbers, i.e., for each inequivalent permutation of the external points, it was 
shown that one can write 

3 

I{u, v) = J2 log' ^ fk{x) + 0{u) , (3.1) 

fc=0 

where I{u,v) denotes either the Easy or the Hard integral, and v = 1 — x + 0{x). The 
coefficients fk{x) were expressed as combinations of terms of the form 

f:^-^S,is) or f:jf^, SAs), (3.2) 

s=l s=l 

where Sj{s) are nested harmonic sums [12], 

S.{s) = ±l- and %(.) = E^- (3-3) 

71=1 n=l 

To arrive at such explicit results for the coefficients fk{x) a kind of experimental math- 
ematics suggested in ref. [31] was applied: the evaluation of the first terms in the expansion 
in X gave a hint about the possible dependence of the coefficient at the n-th power of x. 
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Then an ansatz in the form of a hnear combination of nested sums was constructed and 
the coefficients in this ansatz were fixed by the information about the first terms. Finally, 
the validity of the ansatz was confirmed using information about the next terms. The 
complete x-expansion was thus inferred from the leading terms. 

For the purpose of this paper, it is more convenient to work with polylogarithmic 
functions in x rather than harmonic sums. Indeed, sums of the type fl3.2p can easily be 
performed in terms of harmonic polylogarithms using the algorithms described in ref. |43] . 
We note, however, that during the summation process, sums of the type (13. 2 p with i = 
are generated. Sums of this type are strictly speaking not covered by the algorithms of 
ref. |13], but we can easily reduce them to the case i using the following procedure, 

oo oo s ^ ,oos, 

E S^^i') =-J2^'Y.^S^^) = -T.^'J2-. SAn) , (3.4) 

s=l s=l n=l ^ s=0 n=0 

where the last step follows from 5^(0) = 0. Reshuffling the sum by letting s = rii + n, we 
obtain the following relation which is a special case of eq. (96) in ref. |34j : 

E %(^) = - E E 77 ^» = E ^ s^'^) ■ (3-5) 

s=l ni=0 n=0 s=l 

The last sum is now again of the type (13. 2 p and can be dealt with using the algorithms of 
ref. Hg. 

Performing all the sums that appear in the results of ref. [23], we find for example 

2 2 ^ / 

^13 ^24 = ^ (^-^2,2,1 " -^^2,1,2 + -^^1,3,1 + 2iJl,2,l,l — -f^l,l,3 " 2iJl, 1,1,2 (3.6) 

- 6(3^^2 - 6C3^1,l) - - (2C3i^2,l - 4(3^^1,2 + 4(3 ^^1,1,1 + ^^3,2,1 
~^ -^3,1,2 + -^2,3,1 ~ -^2,1,3 + 2i?i,4,i + 2 i/1,3,1,1 + 2i7i,2,2,l ~ 2iifi,i,4 

- 2i7i,i,2,2 - 2i7i,i,i,3 - 6 C3i^3) + Oiu) , 

^13 ^24 -^12;34 = — (^-^1,1,2,1 — -^^1,1,1,2 " 6^3-^^1,1 j — (^4i72,l,2,l " 4 -^2,1,1,2 (3-7) 

+ 4ifi_i,3,i — -ffi, 1,2,1,1 — 4ifi, 1,1,3 + -^1,1,1,2,1 — 24(3if2,i + 6(3-^^1,1,1 j 

where we used the compressed notation, e.g., -f^2,i,i,2 = H{0,1,1,1,0,1; x). The results 
for the other orientations are rather lengthy, so we do not show them here, but we collect 
them in Appendix [Al Let us however comment about the structure of the functions fk{x) 
that appear in the expansions. The functions /fc(x) can always be written in the form 

h{x) = X PPLs in x] , (3.8) 

I 
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where Rk^i{x) may represent any of the following rational functions 

1 1 1 



X xil — x) 



(3.9) 



We note that the last rational function only enters the asymptotic expansion of -ffi3;24- 

The aim of this paper is to compute the Easy and Hard integrals by writing for each 
integral an ansatz of the form 

Ri{x,x)Pi{x,x), (3.10) 



and to fix the coefficients that appear in the ansatz by matching the limit x — > to 
the asymptotic expansions presented in this section. In the previous section we argued 
that a natural space of functions for the polylogarithmic part Pi{x,x) are functions that 
are single-valued in the complex x plane in Euclidean space. We however still need to 
determine the rational prefactors Ri{x, x), which are not constrained by single- valuedness. 

A natural ansatz would consist in using the same rational prefactors as those appearing 
in the ladder type integrals. For ladder type integrals we have 

i?^'^^-(x,x) = -^— , a en, (3.11) 
[x — XJ" 

plus all possible transformations of this function obtained from the action of the 5*3 sym- 
metry fl2.25p . Then in the limit m — )■ we obtain 

]imR}t'''-{x,x) = ^. (3.12) 

We see that the rational prefactors that appear in the ladder-type integrals can only give 
rise to rational prefactors in the asymptotic expansions with are pure powers of x, and so 
they can never account for the rational function 1 / (x(l— x)) that appears in the asymptotic 
expansion of -ffi3;24- We thus need to consider more general prefactors than those appearing 
in the ladder-type integrals. This issue will be addressed in the next sections. 

4 The Easy integral 

4.1 Residues of the Easy integral 

The Easy integral is defined as 

_ ^23-^24 f ci^xs d'^XQ d^Xj x\q ,^ 

'TT-S / ( /7»2 ^2 /Y>2 ^ /y>2 / /Y>2 ^2 /Y>2 \^2 / ,-y>2 ^2 ,-y>2 \ ^ ' 

n J l,X]^5X25X35jX5g|^X2eX3gX4gjXQ7l^X]^7X27X47j 

To find all its leading singularities we order the integrations as follows 

2 

XoqXi 



7-1 _ ■-'23-^24 
-C^12;34 — ^ 

7r° 



d^Xg xfg f f (i^Xs \ f f d^XY 



^ ^ ^ , , ^ ^ ^ ^ , , , 2222 

rift^ rr>'-' /y^ \ / rr>^ rf-' /y^ I \ I ry-^ rr>^ 

■^26 "''36 -^46 V-' "''15"''25"''35"''56/ \~' "''17"''27"''47"''67 



(4.2) 
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First the Xj and X5 integrations: they are both the same as the massive box computed 
in the Introduction and thus give leading singularities (see eq. f ll.l4p ) 



4A1236 A 1246 



respectively. So we can move directly to the final integration 

^'•^ (4.4) 



(i'^xe x\ 

16 vr^./ x|6a;|6x|6Ai236Al246' 



Here there are five factors in the denominator and we want to take the residues when 
four of them vanish to compute the leading singularity, so there are various choices to 
consider. The simplest option is to cut the three propagators 1/x^g. Then on this cut we 
have Ai236|cut = '^^i^^'iz Ai246|cut = ='=2;^gX24, where the vertical line indicates the value 
on the cut, and the integral reduces to the massive box. This simplification of the A factors 
is similar to the phenomenon of composite leading singularities [H]. Thus cutting either 
of the two As will result iE0 

leading singularity #1 of Ei2;34 = ± „, • (4.5) 

64 7r^Ai234 

The only other possibility is cutting both A's. There are then three possibilities, firstly 
we could cut x\q and x|g as well as the two A's. On this cut A1236 reduces to ix^gXga and 
one obtains residue #1 again. Similarly in the second case where we cut x^g, x\q and the 
two As. 

So finally we consider the case where we cut x\q^ x\q and the two A's. In this case 
Ai236|cut = ±(xfgX23 - xl^xlo) and Ai246|cut = ±{xIqxI^ - xl^xlo). Notice that setting 

„2 



Ai236 = Ai246 = means setting xfg = X2g = 0. We then need to compute the Jacobian 

^36' ^46' A 1236 5 Ai246 



associated with cutting a;|g, x\q, A1236, Ai 



det 



^(^36; ^46^ Al236; Ai246' 
(9Xg 



cut 



ZCiUUei IXgg, X4g, X]^gX23 X23X2g, J']^gX24 ■i^\A,-i^2io 

±16 det (a^36) ^46' ^16' ^26)(^23^14 ^24^13) | cut 
±4Ai234(2;23^14 ~" ^24^13) ' 



(4.6) 

cut ^ ' 



The result of the x^ integral (14.41) is 

1 X 



16 



64 7r6 x\\x2ZA^{AzA^- AaAz) 



(4.7) 



cut 



"'With a slight abuse of language, in the following we use the word 'cut' to designate that we look at 
the zeroes of a certain denominator factor. 
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X 1 oX 



At this point there is a subtlety, since on the cut we have simultaneously x^gXga — 

26 — and so is undefined. More specifically, 



13-^26 



16-^24 



14-^26 



0, i.e. xIq 



X n 



the integral depends on whether we take a^igX^g— x^gXgg = first or x^qxI^—xI^x^q = first. 
So we get two possibilities (after multiplying by the external factors x^^x^^^ in eq. 

leading singularity #2 of -Ei2;34 '- 

leading singularity #3 of -Ei2;34 '- 



™2 fy»2 

. •^13"''24 

"64 vr6 Ai234(a;i3X?4 



/y»2 ^2 

X24X]^3j 



"64 tt'^ Ai234(a;i3X 
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/>! 2 ,-yi2 ^ 

"''24"''13J 



(4.8) 
(4.9) 



We conclude that the Easy integral takes the 'leading singularity times pure function' 
formll 



E 



12;34 



rp 2 /y>2 

"''13"''24 L 



X — x 



V E^^Hx, x) 
_ \ _ ^ ' ^ 

(a; — x){v — 1) (x — x)(f — 1) 



(4.10) 



We note that the 0:3 -H- 0:4 symmetry relates E'^^^ and i?^'^^ Furthermore, putting everything 
over a common denominator it is easy to see that i?^") can be absorbed into the other two 
functions. We conclude that there is in fact only one independent function, and the Easy 
integral can be written in terms of a single pure function E{x,x) as 



E, 



1 



E{x, x) + V E 



X 



X 



a2;34 — 2~r ^V? T\ 

XI3X24 (x — Xj(f — 1) [ \X 

The function E{x, x) is antisymmetric under the interchange of x, x 

E{x, x) = —E{x, x) , 



1' X - 1 



(4.11) 



(4.12) 



to ensure that £^i2;34 is a symmetric function of x, x, but it possesses no other symmetry. 

The other two orientations of the Easy integral are then found by permuting various 
points and are given by 



E 



13;24 



E 



14;23 



^13^24 — x){u — v) 
1 



uE(-,^ 

\X X 

Eil-xA 



+ vE 



1 



1 - x' 1 



X 



X 



uE 1 



X 



(4.13) 
(4.14) 



It is thus enough to have an expression for E{x, x) to determine all possible orientations 
of the Easy integral. The functional form of E{x, x) will be the purpose of the rest of this 
section. 



similar form of the Easy leading singularities, as well as those of the Hard integral discussed in the 
next section, was independently obtained by S. Caron-Huot [45] . 
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4.2 The symbol of E( 

In this subsection we determine the symbol of E{x,x), and in the next section we de- 
scribe its uphft to a function. This strategy seems over-comphcated in the case at hand, 
because E{x,x) can in fact directly be obtained in terms of SVHPLs of weight six from 
its asymptotic expansion using the method described in Section 12.31 The two-step deriva- 
tion (symbol and subsequent uplift) is included mainly for pedagogical purposes because 
it equally applies to the Hard integral and our four-loop example, where the functions are 
not writeable in terms of SVHPLs only so that a direct method yet has to be found. 

Returning to the Easy integral, we start by writing down the most general tensor of 
rank six that 

• has all its entries drawn from the set {x, 1 — x,x,l — x}, 

• satisfies the first entry condition, i.e. the first factors in each tensor are either xx or 

(l-x)(l-a;), 

• is odd under an exchange of x and x. 

This results in a tensor that depends on 2 ■ 4^/2 = 1024 free coefficients (which we assume 
to be rational numbers). Imposing the integrability condition (12. 9p reduces the number of 
free coefficients to 28, which is the number of SVHPLs of weight six that are odd under an 
exchange of x and x. The remaining free coefficients can be fixed by matching to the limit 
ti — >■ 0, — !■ 1, or equivalently x — )■ 0. 

In order to take the limit, we drop every term in the symbol containing an entry 1 — x 
and we replace x — )■ u/x, upon which the singularity is hidden in m. As a result, every 
permutation of our ansatz yields a symbol composed of the three letters {u, x, 1 —x}. This 
tensor can immediately be matched to the symbol of the asymptotic expansion of the Easy 
integral discussed in Section [31 Explicitly, the limits 



^13^24 -^12;34 



lim E[x, x) + lim E \ , _ ''^ 

x-s-O a;->0 V a; — 1 X 



+-limE f^— ^ (4.15) 

X x^O \X - I X - I J 

xls4,Ei3,2, ^ --hmEf-^,-^) (4.16) 

X x^O \1 — X 1 — X J 

x1.xLEi4.23 - lim E(l - X, 1 - x) (4.17) 

X x^O 

can be matched with the asymptotic expansions recast as HPLs. All three conditions are 
consistent with our ansatz; each of them on its own suffices to determine all remaining 
constants. The resulting symbol is a linear combination of 1024 tensors with entries drawn 
from the set {x, 1 — x,x,l — x} and with coefficients {±1, ±2}. 

Note that the uniqueness of the uplift procedure for SVHPLs given in Section [273] implies 
that each asymptotic limit is sufficient to fix the symbol. 
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4.3 The analytic result for E{x,x): uplifting from the symbol 

In this section we determine the function E{x,x) defined in eq. fl4.1ip starting from its 
symboL As the symbol has all its entries drawn from the set {x,l — x,x,l — x}, the 
function E{x,x) can be expressed in terms of the SVHPLs classified in [27] . Additional 
single-valued term^ proportional to zeta values can be fixed by again appealing to the 
asymptotic expansion of the integral. 

We start by writing down an ansatz for E{ linear combination of weight six of 

SVHPLs that is odd under exchange of x and x. Note that we have some freedom w.r.t. 
the basis for our ansatz. In the following we choose basis elements containing a single 
factor of the form Lg_{x). This ensures that all the terms are linearly independent. 

Next we fix the free coefficients in our ansatz by requiring its symbol to agree with that 
of E{x,x) determined in the previous section. As we had started from SVHPLs with the 
correct symmetries and weight, all coefficients are fixed in a unique way. We arrive at the 
following expression for E{x,x): 



For clarity, we suppressed the argument of the L functions and we employed the com- 
pressed notation for HPLs, e.g., 1/3,2,1 = /^o,o,i,o,i,i(2^5 a;). The asymptotic limits of the last 
expression correctly reproduce the terms proportional to zeta values in eq. (13. 7p and the 
formulae in Appendix |X1 

4.4 The analytic result for E{x,x): the direct approach 

Here we quickly give the direct method for obtaining E{x, x) explicitly from its asymptotics 
via the method outlined in Section 12.31 

The asymptotic value of the Easy integral in the permutation -E'i2;34 is given in Ap- 
pendix |X1 Comparing eq. (IA.1|) with eq. fl4.15p and further writing log u = log x + log x 
and expanding out products of functions we find for the asymptotic value of E{x, x): 



where P is a polynomial in logx with coefficients that are HPLs in x. From the discussion 
in Section 12.31 we know that there is a unique combination of SVHPLs with this precise 

^In principle we cannot exclude at this stage more complicated functions of weight less than six multi- 
plied by zeta values. 




2-^3,1,2 + 4L3^2,0 



(4.18) 




(4.19) 
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asymptotic behavior, and so we find a natural ansatz for E{x,x), 

E{x, x) = 4(^3^2, 1 + 2£2,4 — 2/^4 2 + -^1^2,3 ~ -^1,3,2 — '^^1,4,0 + -^2,1,3 ~ '^3,1,2 + 2£3^2,0 
~ '^1,3,1,0 + -^2, 1,2,0 ~ 2^2,2,0,0 " ^2,2,1,0 + -^3,1, 1,0 + 2£l, 2,0,0,0 + -^1, 2, 1,0,0 

- -^2,1,1,0,0 - 20C5/:i + 8C3A + 2(3/^1,2 . (4.20) 

We have hfted this function from its asymptotics in just one hmit x — while we also 
know two other limits of this function given in eq. (13. 7p and Appendix A. Remarkably, 
eq. (l4.2Up is automatically consistent with these two limits, giving a strong indication that 
it is indeed the right function. Furthermore, eq. (14.201) can then in turn be rewritten in 
a way that makes the antisymmetry under exchange of x and x manifest, and we recover 
eq. (14.201) . Note also that antisymmetry in x <H- x was not input anywhere, and the fact 
that the resulting function is indeed antisymmetric is a non-trivial consistency check. 

As an aside we also note here that the form of E{x, x), expressed in the particular basis 
of SVHPLs we chose to work with, is very simple, having only coefficients ±1 or ±2 for 
the polylogarithms of weight six. Indeed other orientations of E have even simpler forms, 
for instance 

1/x) = £2,4 ~ -^3,3 ^ -^1,2,3 + -^1,3,2 — -^1,4,0 ~ -^2,1,3 + -^3,1,2 ~ -^4,0,0 + -^4,1,0 

+ -^1,3,0,0 + -^1,3,1,0 — -^2, 1,2,0 + -^2,2, 1,0 + -^^3,0,0,0 — -^3,1, 1,0 — -^^1,2, 1,0,0 (4.21) 
— -£2, 1,0,0,0 + -^^2, 1,1,0,0 + 8C3'^3 — 2(^3£i^2 — QC'i^2fi — 4C3'C2,1 , 

with all coefficients of the weight six SVHPLs being ±1, or in the manifestly antisymmetric 
form with all weight six SVHPLs with coefficient +1 

E{l/x, 1/x) = ^2,4 + -^^1,3,2 + -^3,1,2 + -^4,1,0 + -^1,3,0,0 (4 22) 

+ ^1,3 ,1,0 + -^2,2,1,0 + -^3,0,0,0 + -^2,1,1,0,0 + 6C3-^3 ~ '^(3^2,1- 



4.5 Numerical consistency tests for E 

We have determined the analytic result for the Easy integral relying on the knowledge of 
its residues, symbol and asymptotic expansions. In order to check the correctness of the 
result, we evaluated -E'i4;23 numericalljlj and compared it to a direct numerical evaluation 
of the coordinate space integral using FIESTA [IHIHH]. 

To be specific, we evaluate the conformally- invariant function 3:^3X24 -Ei4;23. Applying 
a conformal transformation to send X4 to infinity, the integral takes the simplified form, 

lim x^ ^14-23 = — [ d^x^d'^xed'^xj xl^xj^ ^ 

a:4-i'00 71^ J (X]^5X25)x5g (x2gX3g)Xgy (x^^^Xg^) 

with only 8 propagators. We use the remaining freedom to fix x^3 = 1 so that u = x^2 
and V = X23. Other numerical values for x^3 are possible, of course, but we found that this 
choice yields relatively stable numerics. 



^AU polylogarithms appearing in this paper have been evaluated numerically using the GiNaC 
and HPL [37] packages. 
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Analytic 


r ItLO i A 





n 1 
U.i 


n o 


oA.oOOA 


OZ.OOOO 


6.6e-7 


0.2 


0.3 


57.0467 


57.0468 


3.2e-8 


0.3 


0.1 


90.3540 


90.3539 


5.9e-8 


0.4 


0.5 


37.1108 


37.1108 


1.9e-8 


0.5 


0.6 


31.9626 


31.9626 


1.9e-8 


0.6 


0.2 


54.2881 


54.2881 


6.9e-8 


0.7 


0.3 


42.6519 


42.6519 


4.4e-8 


0.8 


0.9 


23.0199 


23.0199 


1.7e-8 


0.9 


0.5 


30.8195 


30.8195 


2.4e-8 



Table 1: Numerical comparison of the analytic result for 0:^3X24 -E'i4;23 against FIESTA for 
several values of the conformal cross ratios. 



After Feynman parameterisation, the integral is only seven-dimensional and can be 
evaluated with off-the-shelf software. We generate the integrand with FIESTA and perform 
the numerical integration with a stand-alone version of CIntegrate. Using the algorithm 
Divonn^, we obtain roughly five digits of precision after five million function evaluations. 

In total, we checked 40 different pairs of values for the cross ratios and we found very 
good agreement in all cases. A sample of the numerical checks is shown in Table [1] Note 
that 6 denotes the relative error between the analytic result and the number obtained by 
FIESTA, 

Nanalytic — NpiESTA 



analytic 



+ NpiESTA 



(4.24) 



5 The Hard integral 

5.1 Residues of the Hard integral 

To find all the leading singularities we consider each integration sequentially as follows 



X 



TT _ -34 
n^o■?,A — — — 



2222 

"''16"''26-^36"''46 L 



Let US start with the x-j integration, 



(i^Xs x\ 



1111 

ry^ ry ^ ry^ ry ^ 

"''15"'' 25-^35 "''45 



d^Xj 



2 2 2 2 

ry ^ ry^ ry*-' ry*-' 



(5.1) 



d'^xj 



1111' 

ry'-' ry'-' ry'-' ry'-' 

X2yX^Y-^57-^67 



(5.2) 



^Experience shows that Divonne outperforms other algorithms of the Cuba hbrary for problems roughly 
this size. 
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This is simply the off-shell box considered in Section [H and so its leading singularities are 
(see eq. firTi]) ) 

± • (5.3) 

4 A3456 



Next we turn to the integration, which now takes the form 

(5.4) 



d X5 Xgg 

^15^25^35^45-^3456 

There are five factors in the denominator, and we want to cut four of them to compute 
the leading singularity. The simplest option is to cut the four propagators l/a;^5. Doing 
so would yield a new Jacobian factor I/A1234 (exactly as in the previous subsection) and 
freeze A3456|cut = '^^le^'zi- This latter factor simply cancels the numerator and we are left 
with the final Xq integration being that of the box in the Introduction. Putting everything 
together, the leading singularity for this choice is 

leading singularity #1 of -f/'i2;34 = ± „. — • (5-5) 

o47r A1234 

Returning to the X5 integration, eq. (15. 4p . we must consider the possibility of cutting 
A3456 and three other propagators. Cutting x\^ and x\^ immediately freezes A3456|cut = 
±xIqxI^ which is canceled by the numerator. Thus it is not possible to cut these two 
propagators and A3456. However, cutting a^is, a^25' ^35 -^3456 is possible (the only other 
possibility, i.e. cutting xl^,X2^,xl^ and A3456, gives the same result by by invariance of the 
integral under exchange of X3 and X4). Indeed one finds that when = 0, 

A3456 = ±(a;45a;36 - a;56^34) • (5-6) 

To compute the leading singularity associated with this pole we need to compute the 
Jacobian 

^(^15' ^25' ^35' -^3456) 



J ^ ^ ^^^ ^^^7^^ , (5.7) 

As in the box case, it is useful to consider the square of J (on the cut), 

J - 16det 1^ _ ^^^^^^^^^^ (9A3456/9x5)2 ) ' ^^'^^ 
The result of the X5 integration is then simply 

X 



^56 



7^2 



^36 



7^2 
cut "'•^34 



cut 



(5.9) 



where the second equality follows since x^q and X45 are to be evaluated on the cut (indicated 
by the vertical line) for which xI^x'^q — x^qX^^ = 0. Finally we need to turn to the remaining 
xq integral. We are simply left with 



1 f d'^XQ 



IQn^ J xIqX2qxIqJ 



(5.10) 

cut 
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where we note that the x|g propagator term has canceled with the numerator in eq. fl5.9l) . 
So we have no choice left for the quadruple cut as there are only four poles. In fact on the 
other cut of the three propagators we find J|cut = '^{xi4,X23 ^ ^13^24)^36' ^o this brings 
back the propagator x"^ 



36- 



Computing the Jacobian associated with this final integration thus yields the final result 
for the leading singularity, 



leading singularity #2 of Hi2- 



1 



34 



"647r6 {xj^xl^ - xl3xl^)Xi234 ' 



(5.11) 



We conclude that the Hard integral can be written as these leading singularities times 
pure functions, i.e. it has the form 



H 



12;34 



4 4 
■^13-^24: L 



X, X 



XI 



(x — x)2 {v — l){x — x) 



(5.12) 



where are pure polylogarithmic functions. The pure functions must furthermore 

satisfy the following properties 



i7(")(x,x) = i7(")(x,x), 

ff(«)(a;,x) = H^^\x/{x - l),x/(x - 1)) 



X, X 



(5.13) 



H^''\x,x) = H^''\x/{x-l),x/{x-l)) 



in order that i?i2;34 be symmetric in x, x and under the permutation Xi 0x2. Furthermore 
we would expect that H^""^ (x, x) = in order to cancel the pole at x — x. In fact it will 
turn out in this section that even without imposing this condition by hand we will arrive 
at a unique result which nevertheless has this particular property. 
By swapping the points around we automatically get 



H 



13;24 



H 



14:23 



4 4 
•^13-^24 



^13^24 



i/('^)(l/x, 1/x) ^ H^''\l/x,l/x) 



X — X 



H^°-Xl -X, 1 



X 



(m — w)(x — X 
/7W(1 



X, 1 



x) 



X 



x)^ 



m)(x — x) 



(5.14) 
(5.15) 



5.2 The symbols of H^^^x^x) and H^^\x,x) 

In order to determine the pure functions contributing to the Hard integral, we proceed just 
like for the Easy integral and first determine the symbol. For the Hard integral we have 
to start from two ansatze for the symbols S[H^°'\x^ x)] and S[H^^\x^ x)]. While both pure 
functions are invariant under the exchange Xi o X2, S[H^'^'>] must be symmetric under the 
exchange of x, x and S[H^^''] has to be antisymmetric, cf. eq. (15.131) . Going through exactly 
the same steps as for E we find that the single-variable limits of the symbols cannot be 
matched against the data from the asymptotic expansions using only entries from the set 
{x, 1 — X, X, 1 — x}. We thus need to enlarge the ansatz. 

Previously, the letter x — x ~ A1234 has been encountered in ref. [301150] in a similar 
context. We therefore consider all possible integrable symbols made from the letters {x, 1 — 
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x,x,l — x,x — x} which obey the initial entry condition fl2.13p . In the case of the Easy 
integral, the integrability condition only implied that terms depending on both x and x 
come from products of single-variable functions. Here, on the other hand, the condition is 
more non-trivial since, for example, 

X 

d log — Ad log(x — x) = d log x Ad log x , 

1-x (5-16) 
d log A d log(x — x) = d log(l — x) A d log(l — x) . 

1 — X 

We summarise the dimensions of the spaces of such symbols, split according to parity under 
exchange of x and x, in Table [H 



Weight Even Odd 

1 2 

2 3 1 

3 6 3 

4 12 9 

5 28 24 

6 69 65 

Table 2: Dimensions of the spaces of integrable symbols with entries drawn from the set 
{x, 1— x,x,l — x,x — x} and split according to the parity under exchange of x and x. 

Given our ansatz for the symbols of the functions we are looking for, we then match 
against the twist two asymptotics as described previously. We find a unique solution for 
the symbols of both H^°'^ and H^''^ compatible with all asymptotic limits. Interestingly, 
the limit of -ffi3;24 leaves one undetermined parameter in SlH^""^], which we may fix by 
appealing to another limit. In the resulting symbols, the letter x — x occurs only in the 
last two entries of SlH^""^] while it is absent from S[H^''^]. Although we did not impose 
this as a constraint, S[H^'^^] goes to zero when x x, which is necessary since the integral 
cannot have a pole ai x = x. 

5.3 The analytic results for H^°'\x^x) and H^^\x,x) 

In this section we integrate the symbol of the Hard integral to a function, i.e. we determine 
the full answers for the functions H^"'\x,x) and H^'^\x,x) that contribute to the Hard 
integral Hi2-u- 

In the previous section we already argued that the symbol of H^''\x, x) has all its entries 
drawn form the set {x, 1 — x, x,l — x}, and so it is reasonable to assume that H^''\x,x) 
can be expressed in terms of SVHPLs only. We may therefore proceed by lifting directly 
from the asymptotic form as we did in Section 14.41 for the Easy integral. By comparing the 
form of -ffi3;24, eq. (15.141) . with its asymptotic value fll.l4p we can read off the asymptotic 
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form of H{l/x, 1/x). Writing logM as logx + logx, expanding out all the functions and 
neglecting logx terms, we can the lift directly to the full function by simply converting 
HPLs to SVHPLs. In this way we arrive at 

H'^'^l/x, 1/x) = 2£2,4 - 2£3,3 - 2£i,i,4 - 2£i,4,o + 2£i,4,i - 2£2,3,i + 2£3,i,2 
~ 2/^4^0,0 + 2£4,i,o + 2£i^i^i^3 + 2£i 1 3 + 2£i 3 o,o ~ 2£i 3 1 1 
~ 2£2, 1,1,2 + 2/^2,1,2,1 + 2/^3^0,0,0 — 2/23^1^1^0 — 2/2i, 1^1^2,1 — 2/2i^i,2,i,o (5-17) 
+ 2/2i,i,2,i,i — 2£i^2, 1,0,0 + 2/2i,2,i,i,o ^ 2/22,1,0,0,0 + 2/22, 1,1,0,0 

+ I6C3/23 - 16C3>C2,1 . 

Other orientations although still quite simple do not all share the property that they only 
have coefficients ±2. Using the basis of SVHPLs that makes the parity under exchange of 
X and X explicit, we can write the last equation in the equivalent form 

H^''\X, X) = 16^2,4 - 16L4,2 - 8Li,3,2 - 8Li,4,l + 8L2,l,3 - 8^2,2,2 + 8^2,3,1 - 8^3,1,2 
+ 16/^3^2,0 + 8-^3^2,1 — 8L4^i^i + 4Li^2,2,l — 8Li^3^i^i — 41/2,1,1,2 + 8-^2,1,2,1 

— 8L2,2,i,o — 4L2,2,i,i + 8/73^1^1^0 — 4Li_i^2,i,i — 24L2,i,i,i,o • (5.18) 

Next, we turn to the function H^"'\x,x). As the symbol of H^"'\x,x) contains the 
entry x — x, it cannot be expressed through SVHPLs only. Single-valued functions whose 
symbols have entries drawn form the set {x, 1— x,x,l — x,x — x} have been studied 
up to weight four in ref. [30], and a basis for the corresponding space of functions was 
constructed. The resulting single- valued functions are combinations of logarithms of x and 
X and multiple polylogarithms G(ai, . . . , a„; 1), with G {0, Note that the 
harmonic polylogarithms form a subalgebra of this class of functions, because we have, 
e.g., 

g(o,-,--i) =//(0,1,1;x). (5.19) 
\ X X / 

This class of single-valued functions thus provides a natural extension of the SVHPLs we 
have encountered so far. In the following we show how we can integrate the symbol of 
H^'^\x,x) in terms of these functions. The basic idea is the same as for the case of the 
SVHPLs: we would like to write down the most general linear combination of multiple 
polylogarithms of this type and fix their coefficients by matching to the symbol and the 
asymptotic expansion of H^^\x^x). Unlike the SVHPL case, however, some of the steps 
are technically more involved, and we therefore discuss these points in detail. 

Let us denote by Q the algebra generated by logx and logx and by multiple polylog- 
arithms G(ai, . . . , ttri'i 1), with tti G {0, 1/x, 1/x}, with coefficients that are polynomials in 
multiple zeta values. Note that without loss of generality we may assume that a„ 7^ 0. 
In the following we denote by Q"^ the linear subspaces of Q of the functions that are re- 
spectively even and odd under an exchange of x and x. Our first goal will be to construct 
a basis for the algebra as well as for its even and odd subspaces. As we know the 
generators of the algebra we automatically know a basis for the underlying vector space 
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for every weight. It is however often desirable to choose a basis that "recycles" as much as 
possible information from lower weights, i.e. we would like to choose a basis that explicitly 
includes all possible products of lower weight basis elements. Such a basis can always eas- 
ily be constructed: indeed, a theorem by Radford [51] states that every shuffle algebra is 
isomorphic to the polynomial algebra constructed out of its Lyndon words. In our case, we 
immediately obtain a basis for Q by taking products of log x and \ogx and G{ai, . . . , a„; 1), 
where (ai, . . . , a„) is a Lyndon word in the three letters {0, Next, we can eas- 
ily construct a basis for the eigenspaces Q"^ by decomposing each (indecomposable) basis 
function into its even and odd parts. In the following we use the shorthands 

Gt,,...,^S^u ■ ■ ■ , Xfc) = \g[^^_^. ^, . . . , 0^_^, ^; l) ± (x O x) . (5.20) 

mi— 1 mfc— 1 

In doing so we have seemingly doubled the number of basis functions, and so not all the 
eigenf unctions corresponding to Lyndon words can be independent. Indeed, we have for 
example 

GUx,x) = \Gt{xf-\G~,{xf. (5.21) 

It is easy to check this relation by computing the symbol of both sides of the equation. 
Similar relations can be obtained without much effort for higher weight functions. The 
resulting linearly independent set of functions are the desired bases for the eigenspaces. 
We can now immediately write down the most general linear combination of elements of 
weight six in and determine the coefflcients by matching to the symbol of H^'^\x,x). 
As we are working with a basis, all the coefficients are fixed uniquely. 

At this stage we have determined a function in whose symbol matches the symbol of 
H^"'\x, x). We have however not yet fixed the terms proportional to zeta values. We start 
by parametrising these terms by writing down all possible products of zeta values and basis 
functions in Q~^. Some of the free parameters can immediately be fixed by requiring the 
function to vanish for x = x and by matching to the asymptotic expansion. Note that our 
basis makes it particularly easy to compute the leading term in the limit x — 0, because 

limG±(...,x,...)=0. (5.22) 

In other words, the small u limit can easily be approached by dropping all terms which 
involve (non-trivial) basis functions that depend on x. The remaining terms only depend 
on logx and harmonic polylogarithms in x. However, unlike for SVHPLs, matching to 
the asymptotic expansions does not fix uniquely the terms proportional to zeta values. 
The reason for this is that, while in the SVHPL case we could rely on our knowledge of 
a basis for the single-valued subspace of harmonic polylogarithms, in the present case we 
have been working with a basis for the full space, and so the function we obtain might still 
contain non-trivial discontinuities. In the remainder of this section we discuss how on can 
fix this ambiguity. 

In ref. [30] a criterion was given that allows one to determine whether a given func- 
tion is single-valued. In order to understand the criterion, let us consider the algebra 
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Q generated by multiple polylogarithms G(ai, . . . , a„; a„+i), with G {0, 1/x} and 
a„_|_i G {0, 1, with coefficients that are polynomials in multiple zeta values. Note 
that Q contains ^ as a subalgebra. The reason to consider the larger algebra Q is that Q car- 
ries a Hopf algebra structurqj [52], i-e. G can be equipped with a coproduct A : Q ^ Q ^Q. 
Consider now the subspace Qsv oi Q consisting of single-valued functions. It is easy to see 
that Qsv is a subalgebra of Q. However, it is not a sub-Hopf algebra, but rather Qsv is a 
^-comodule, i.e. A : Qsv ^ Qsv ® Q- In other words, when acting with the coproduct on 
a single-valued function, the ffist factor in the coproduct must itself be single-valued. As 
a simple example, we have 

A(L2) = -Lo®log— ^ + -Li®log-. (5.23) 
2 1 — X 2 X 

Note that this is a natural extension of the ffist entry condition discussed in Section [2l 
This criterion can now be used to recursively fix the remaining ambiguities to obtain a 
single-valued function. In particular, in ref. [SU] an explicit basis up to weight four was 
constructed for Qsv- We extended this construction and obtained a complete basis at 
weight five, and we refer to ref. [30] about the construction of the basis. All the remaining 
ambiguities can then easily be fixed by requiring that after acting with the coproduct, the 
ffist factor can be decomposed into the basis of Qsv up to weight five. We then finally 
arrive at 

28 136 160 

H^-){X,X) = H{X,X) - yC3^1,2 + 164(3^2,0 + —(3^2,1 " —^^3^2,1 " QQLqL.^A 

- ^^0^2,3 + ^-^2^3,1 + yLoL3,2 + 16LiL3,2 + 36LoL4,i + 64LiL4,i 

TO 26 
+ -^LqLi 2^2 + '24:LqLi '3 i + —LiLi ;i i — 8/72/72,1,1 + 64/79/72,1,2 

58 50 88 

—LoL2,2,0 — 4/79/72,2,1 + — /^l/^2,2,l — 12/70/73,1,0 ^-^0-^3,1,1 

32 166 
+ 18/71/73,1,1 — —LqLi i 2,i — 18Lo/7i,2,i,i H — ^/7oL2, 1,1,0 — 8/79/72,1,1,1 

+ 328(3/^3 + 3213^ - UL2U ■ 

The function x) is a single- valued combination of multiple polylogarithms that cannot 
be expressed through SVHPLs alone, 

1-L{x,x) = -128G+2 - 512G'+j - 64Gj^ 2 + 64Gjj 2 - 64^^^ 5 - 128Gj2,i 

+ 64G+^ J - 64G+j^, - 448G+J J + 64^2+^^2,1 + 64Gjj 5,1 + 64G+2,i,i + 64^2+2,1,1 
- 64Gj2,i,i + 128(7^2,1,1 + 128G'j2,i,i + 256Gj^ j + 128G+^ j - 1286*^^ ^^ 



+ 192Gjj^^_^ - 64Gjj^^ J - 64Gjj + 192Gjj j j + I28//2+4 - I28//4+2 
640 „^ 64^+ 256 



+ —Hti,3 - -Hts,! - ^^^1,2 + 64//2ti,i,2 - 64//2+ ,1,1 + 64LoGt2 (5.25) 



^Note that we consider a slightly extended version of the Hopf algebra considered in ref. that allows 
us to include consistently multiple zeta values of even weight, see ref. [551151] . 
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+ 256LoGjj + 32LoGji 2 + 64LoGj2,i + ^^^0^X1,1 + 32i^oGjj ^ + 96LoGjj 

- 32LiGj2,i - SOLiGj^ J - 16LiGjj^^ - 16LiGjj ^ - 64^2^"^ j + 64^4^"^ 

32 

+ 32L2,2G~j - —H^H+^ - QAH+H+^ ^ - 128H^H+ - 64/7f LqG" j j 

- 32L'^G^ ^ — 32Z7qG^j j + 32Z;Q(jr^-|^ 1 1 "I" '^^LiLqG^ ^ + 16ZyxijoC'^2^]^ j 

+ 16LiLoG'2+j J - y ffr^o^2,2 - 48ifrLoL2,i,i + 12H^L,L2,2 + 16L^if2+2 

640 

+ 32X2/^2^1,1 - 64i7f L4L0 + IQH^LiU + 64L3G+1 j - —H^H^^^ 
+ 64(if2";i)^ + 128(iJ3-)2 + 32L0L2G-J - 32LoL2G-^ J - I6L1L2G-J 

H — —LQL2H2 I + 16i/]^ L2iv2,i —H2 LqL2,i — 8H2 LiL2^i — 32H^ LqL2 

- 48iff L3L2 + ?,2H+LqL^ + WH^L^L^ + 32i7f L^G^j - 16iJf LiLqG^j 

+ "^-^o^2ii + "^-^o^tiJ ~ SLiLqG^j - 8LiLoG+^ J + -^-f^i -^0-^2,1 

- 16(i/r)2LoL2,i - 32H^H^Ll + ^(ifr)'^3^o - l2{H^fL^L; + 28/72^^2 

H ^g"^"^ 16LqL2G'j^j — 8Lo-LiL2G'j^j + —H{H2LqL2 — 8H^ H2L1L2 

+ 8{H^fLl + 8(if2+)2L^ + 8{HtfL,L, + |(ifr)2if2^L^ - A{H^fHtL,L, 

160 52 

- 9QH2iH{fLo + — (iJr)'^oi:2 + -^H^LlL2 + AH^ ULIL2 

2 148 
+ AH^ LIL1L2 + H2L0LI + -H2LILI — 8-^2^X0X1 + -^(iff )''Lo 

+ y (i/D^L^ + ^{H^fLlLl - ^{H^YLlL^ - 128(3^^ - 128(3^+,^ 

16 64 
+ yC3(i/-)2Lo + 24C3(//r)'i:i + yC3i^ri^2 , 

where we used the obvious shorthand 

= ^H^ix) ± (x O x) . (5.26) 

and similarly for G^. In addition, for the position of x is indicated by the bars in the 
indices, e.g., 

^t2,3 = ^t2,3(^,^,^)- (5-27) 

Note that we have expressed x) entirely using the basis of constructed at the 
beginning of this section. As a consequence, all the terms are linearly independent and 
there can be no cancellations among different terms. 
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u V Analytic FIESTA 6 



0.1 0.2 269.239 269.236 6.4e-6 

0.2 0.3 136.518 136.518 1.9e-6 

0.3 0.1 204.231 204.230 1.3e-6 

0.4 0.5 61.2506 61.2505 5.0e-7 

0.5 0.6 46.1929 46.1928 3.5e-7 

0.6 0.2 82.7081 82.7080 7.4e-7 

0.7 0.3 57.5219 57.5219 4.7e-7 

0.8 0.9 24.6343 24.6343 2.0e-7 

0.9 0.5 34.1212 34.1212 2.6e-7 



Table 3: Numerical comparison of the analytic result for xf^X24^ -f^i3;24 against FIESTA for 
several values of the conformal cross ratios. 

5.4 Numerical consistency checks for H 

In the previous section we have determined the analytic result for the Hard integral. In 
order to check that our method indeed produced the correct result for the integral, we 
have compared our expression numerically against FIESTA. Specifically, we evaluate the 
conformally-invariant function xf^x H 1^-24 ■ Applying a conformal transformation to send 
X4 to infinity, the integral takes the simplified form. 



with 9 propagators. As we did for -E'i4;23, we use the remaining freedom to fix xf^ = 1 so 
that u = — ^235 perform the numerical evaluation using the same setup. We 

compare at 40 different values, and find excellent agreement in all cases. A small sample 
of the numerical checks is shown in Table |3l 

6 The analytic result for the three-loop correlator 

In the previous sections we computed the Easy and Hard integrals analytically. Using 
eq. fll.8p . we can therefore immediately write down the analytic answer for the three-loop 
correlator of four stress tensor multiplets. We find 



lim i?13;24 = 



/ 



(5.28) 



(^2 Ay«2 ^2 \ ^2 f ^2 "\ ^2 / ^2 fy»2 ^2 "\ ' 
'''15"''25"''35J'^56l"''36J"''67l'''17"''27"''37J 




1 — X 



1 



) 



(6.1) 
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X — X 
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V - 1 
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E(x) 



E 



X 
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u 



Eil-- 

X 



E(l-x) 



u 



u — v 



X 



E 



u — v 



X 



(x — x)" 



(1 + v) H^^Ux) + (1 + n) H^"^ (l-x) + (u + v) H^"^ - 

X 



+ 



X — X 



V - 1 



1 



u 



X) 



X) 



u 



U + V 



u — v 



The pure functions appearing in the correlator are defined in eqs. fl2.3p . fl4.20p . f l5.17p 
and f l5.24p . For clarity, we suppressed the dependence of the pure functions on x, i.e. we 
write f^^\x) = f^^\x,x) and so on. All the pure functions can be expressed in terms 
of SVHPLs, except for H^""^ which contains functions whose symbols involve 
entry. We checked that these contributions do not cancel in the sum over all contributions 
to the correlator. 



7 A four-loop example 

In this section we will discuss a four-loop integral to illustrate how our techniques can be 
applied at higher loops. The example we consider contributes to the four-loop four-point 
function of stress-tensor multiplets in A/" = 4 SYM. Specifically, we consider the Euclidean, 
conformal, four-loop integral. 



- — f d'^x^d'^XQd'^xid'^xsxl^xl^xl^ _ 1 , , , 

14,23 ■ttS / /Y>2 ^2 /Y>2 ^2 /y>2 ^2 ^2 ^2 ^2 ^2 ™2 ^2 /y»2 ^2 ' ' 

n J Xi5X]^8X25X26X37X3gX45X4gX47X48J.5gXg7X78 X]^3X24 

where the cross ratios u and v are defined by eq. (11. 111) . As we will demonstrate in the 
following sections, this integral obeys a second-order differential equation whose solution is 
uniquely specified by imposing single-valued behaviour, similar to the generalised ladders 
considered in ref. [35]. 

The four-loop contribution to the stress-tensor four-point function in A/" = 4 SYM 
contains some integrals that do not obviously obey any such differential equations, and with 
the effort presented here we also wanted to learn to what extent the two-step procedure of 
deriving symbols and subsequently uplifting them to functions can be repeated for those 
cases. Our results are encouraging: the main technical obstable is obtaining sufficient data 
from the asymptotic expansions; we show that this step is indeed feasible, at least for 
and present the results in Section 17.11 Ultimately we find it simpler to evaluate /^^^ by 
solving a differential equation, and in this case the asymptotic expansions provide stringent 
consistency checks. 
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7.1 Asymptotic expansions 

Let us first consider the limits of tlie four-loop integral (17.11) and its point permutations 
for Xi2,xs4 — > 0. We derive expressions for its asymptotic expansion in the limit where 
n — )■ 0, i; — 7- 1 similar to those for the Easy and Hard integrals obtained in Section |3l The 
logarithmic terms can be fully determined, while the non- logarithmic part of the expansion 
requires four-loop IBP techniques that allow us to reach spin 15. This contains enough 
information to fix the C„log'^(n) terms (important for beyond-the-symbol contributions) 
while the purely rational part of the asymptotic series remains partially undetermined. 
However, our experience with Easy and Hard has shown that each of the three coincidence 
limits is (almost) sufficient to pin down the various symbols. Inverting the integrals from 
one orientation to another ties non-logarithmic terms in one expansion to logarithmic ones 
in another, so that we do in fact command over much more data than it superficially seems. 
It is also conceivable to take into account more than the lowest order in u. 

We start by investigating the asymptotic expansion of the integral /i4^23 whose coinci- 
dence limit Xi2,X34 — diverges as log^w. There are three contributing regions: while in 
the first two regions the original integral factors into a product of two two-loop integrals or 
a one-loop integral and a trivial three-loop integral, the third part corresponds to the four- 
loop 'hard' region in which the original integral is simply expanded in the small distances. 
The coefficients of the logarithmically divergent terms in the asymptotic expansion, i.e. 
the coefficients of log^ u and log u, can be worked out from the first two regions alone. It 
is easy to reach high powers in x and we obtain a safe match onto harmonic series of the 
type (13. 2p with i > 1. Similar to the case of the Easy and Hard integrals discussed in 
Section |3l we can sum up the harmonic sums in terms of HPLs. Note that the absence of 
harmonic sums with i = 1 implies the absence of HPLs of the form Hi {x). 

In the hard region, we have explicitly worked out the contribution from spin zero 
through eight, i.e., up to and including terms of 0{x^). By what has been said above about 
the form of the series, this amount of data is sufficient to pin down the terms involving 
zeta values, while we cannot hope to fix the purely rational part where the dimension of 
the ansatz is larger than the number of constraints we can obtain. The linear combination 
displayed below was found from the limit x — of the symbol of the four-loop integral 
derived in subsequent sections. Its expansion around x = reproduces the asymptotic 
expansion of the integral up to 0{x^). We find 

^13 ^24 -^14;23 — ■'^) 

— log^ U -^2,1,3 — -^2,3,1 + -f^3,l,2 — -^^3,2,1 + 2i72, 1,1,2 " 2/^2,2,1,1 + C^i^Hj, + 6if2,l) + 
^ 3j L 

- logM [-4i72,l,4 + 4i/2,4,l - 3i/3,l,3 + 3/^3,3,1 ' 3^^4,1,2 + 3/^4,2,1 ' 4/^2,1,1,3 ' 4/^2,1,2,2 

X L 

+ 4/^2,2,2,1 + 4i/2,3,l,l — 2if3_i^i^2 + 2/^3,2,1,1 + C3( — 18-^^4 — 8/^2,2 " 2/^3, i + 8/^2,1,1 
lOi/2,1,5 + 2-^2,2,4 — 2iJ2,3,3 " 10if2,5,l + 8-^3,1,4 — ^H^^^^i + 6-^4,1,3 — 6iJ4,3,i 



+ 



X 
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+ 6-^5,1,2 — 6-^5,2,1 + 8-^2,1,1,4 + 6/^2, 1,2,3 + 8-^2,1,3,2 — 2-f/'2,l,4,l + 2-^2,2,2,2 " 4-^2,2,3,1 

— 4-ff2,3,l,2 — lO-f/2,3,2,1 — 4-ff2,4,l,l + 4-^3,1,1,3 + 6-^3,1,2,2 " 6-ff3,2,2,l " 4-ff3,3,i,i + 4-^2,1,1,2,2 

— 4-ff2,l,2,2,l — 4-f/'2,2,l,l,2 + 4-^2,2,2,1,1 + Csi'^QHr^ + 8-^/2,3 + 12-ff3,2 — 1211^^1 — 4-^2,1,2 

— 16i72,2,l - 8i/3,l,l) + C5(10f^3 + lOi/2,1)] + Oiu) . 

Next we turn to the asymptotic expansion of the orientation /i2-34- Here the Euchdean 
coincidence hmit xi2,X34 — )■ is finite, and thus the only region we need to analyse is the 
four-loop hard region, for which we have determined the asymptotic expansion up to and 
including terms of 0{x^^). Just like for the non-logarithmic part in the asymptotic expan- 
sion of /j4^23 5 61- (17. 2p . we have fixed the terms proportional to zeta values by matching 
an ansatz in terms of HPLs onto this data, and once again, the terms not containing zeta 
values are taken from the relevant limit of the symbol. We find 

^13 ^24 -^12-34 — (7-3) 



4-^1,3,4 — 4-f/'i,5,2 + 2-^1,1,2,4 — 2-^1,1,4,2 + 2-f/'i,2,l,4 " 2-f/'i,2,3,2 + 2-^1,3,1,3 + 2-^/1,3,3,1 

— 2-ffi,4,i,2 — 2-ffi, 5,1,1 + -f^l, 1,2,1,3 + -f^l, 1,2,3,1 — -f^l,l,3,l,2 " -f^l,l,3,2,l + -f^l,2,l,l,3 + -f^l, 2,1,3,1 
"~ -^^1,2,2,1,2 + -f^l,2,2,2,l — 2-^1,2,3,1,1 + -^1,3,1, 2,1 — -^^1,3,2,1,1 + -f^l, 2, 1,1,2,1 — -^^1,2,1,2,1,1 + 

(3(8^/1,1,3 - 8i/i,2,2 + 4/^1,1,2,1 - 4//i,2,i,i) + 70 CrHi] + Oiu) . 

The expansion around x = of this expression reproduces the asymptotic expansion of the 
integral up to 0{x^^). 

The most complicated integrals appearing in the asymptotic expansion of l['^.23, -^12-34 
are four-loop two-point dimensionally regularised (in position space) integrals which belong 
to the family of integrals contributing to the evaluation of the five-loop contribution to the 
Konishi anomalous dimension [221. 



, f d'^XQdf'-X'jd'^x^d'^x^ 

G(ai, . . . ,014) = j - 



[xi^Y^ [xIjY^ [xi^y-^ [xi^Y^ [xlY- {x^Y' {^lY' 

{^liY' {xl%Y'° i^mY'' i^isY'' (a;?9)"" (a^ig)'^" ' ^'^■^^ 

with various integer indices ai, . . . , ai4 and d = 4 — 2e. 

The complexity of the IBP reduction to master integrals is determined, in a first ap- 
proximation, by the number of positive indices and the maximal deviation from the corner 
point of a sector, which has indices equal to or 1 for non-positive and positive indices, cor- 
respondingly. This deviation can be characterised by the number Yli^u+^'^i ~ '^i^v- 
where i'± are sets of positive (negative) indices. So the most complicated (for an IBP 
reduction) integrals appearing in the contribution of spin s to the asymptotic expansion in 
the short-distance limit have nine positive indices and the deviation from the corner point 
is equal to 2s — 2. It was possible to get results up to spin 15. 
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As in ref. [22] the IBP reduction was performed by the C++ version of the code FIRE |4U] . 
The master integrals of this family either reduce, via a dual transformation, to the cor- 
responding momentum space master integrals [55l[56] or can be taken from ref. [22]. To 
arrive at contributions corresponding to higher spin values, FIRE was combined with a 
recently developed alternative code to solve IBP relations LiteRed [57] based on the alge- 
braic properties of IBP relations revealed in ref. [SB]. (See ref. [SH] where this combination 
was presented within the Mathematica version of FIRE.) 



7.2 A differential equation 

We can use the magic identity [23] on the two-loop ladder subintegral 



/ rfi'^ /7»2 ^2 ,-y>2 ^2 ,-y>2 ^2 
/I J Xi5X25X26"f'45X4gXg7X5g 



The magic identity reads 

I^'^\xi,X2,Xa,Xi) = l'^'^\x2,Xi,X7,X4) , (7.6) 

and using it on the four-loop integral we find 

f d^^rd^Xs (2). ^ ^ ^ ^ 

-'l4;23 -4/^2 ™2 ™2 ™2 ™2 ™2 l-^l, -^2, ^4, J^7j 

_ ]_ f d'^X7d^X8 (2)( N 

I rfi-^ lyt^ rp^ rv*^ lyt^ rv*^ 

li J '^i8'^37'^38*^47'^ 48*^78 

1 /* d x^d x^d X'jd x^x-^jX-^^x^^ . . 

■Tr8 / /y>2 /y>2 ^2 ,-y>2 ^2 ,-y>2 ^2 ,-yi2 ^2 ,-y>2 ^2 ,-y>2 ^2 ^ 

The resulting integral (17. 7p is 'boxable', i.e. we may apply the Laplace operator at the 
point X2. The only propagator which depends on X2 is the one connected to the point x^ 
and we have 



1 



□2— = -47r^54(x25) . (7.^ 

^25 



The effect of the Laplace operator is therefore to reduce the loop order by one [2S]. Thus 
on the full integral /^^^ we have 



U2Jl4;23 - 



d^x^db^Xid^x^XYiX^^x^/^ 



fTr\j I rr>^ nn" ly'^ nn" ly'^ nn" ly'^ rr>" nn^ /y" ry^ /y" 

/I J X;^8X37X38X47X48X78X]^2"''16"''72'''76"''64"''26 

= -4-#^i^l4;23 , (7.9) 
•^12'''24 

where we have recognised the Easy integral, 

_ 1 /" d'^ Xgd^ Xjd'^ xsx'^^X24^XiY ^ ^ f ( \ /"7 1 n^ 

-£^14;23 — — ^ / -2 2 2 2 2 2 2 2 2 2 2~ ~ ^ ^Je[U,V . 

/I J J'i8-^37J'38X47J'48-^78-^16-^72-^76-^64-^26 •^13-^24 
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The differential equation fl7.9p becomes an equation for the function /, 

□2-^/(«, V) = V) . (7.11) 

Applying the chain rule we obtain the following equation in terms of u and f , 

A(^)/(u,i;) = -^/^(u,t;), (7.12) 

where 

A(2) = 4[2(9, + d,) + udl + vdl -{l-u- v)dud,] . (7.13) 
In terms of {x,x) we have 

xxd^dsfix, x) = -/ij(x, x) , (7.14) 

where 

f{x,x) = -{x-x)f{u,v) (7.15) 

and similarly for Je- Note that f{x,x) = —f{x,x). Now we recall that the function 
fE{u,v) defined by eq. fl7.10p in the orientation -Ei4;23 is of the form 

fE{u,v) = - -^^ -\e{1-x,1-x)+xxE(i--,1--)] . (7.16) 

[X — X){1 — XX) L \ X XJ\ 

Hence we find the following equation for /, 

{1 — xx)xxdj,dxf{x,x) = — E{\ — x,\ — x) \ xxE{\ , 1 — — ^ . (7.17) 



X xJ J 

Without examining the equation in great detail we can immediately make the following 
observations about /. 

• The function / is a pure function of weight eight. From eq. fl7.15p the only leading 
singularity of the four-loop integral /^^^ is therefore of the \j{x — x) type, just as for 
the ladders. 

• The final entries of the symbol of /(x, x) can be written as functions only of x or of x, 
but not both together. This follows because the right-hand side of eq. (17.171) contains 
only functions of weight six, whereas there would be a contribution of weight seven 
if the final entries could not be separated into functions of x or x separately. 

• The factor (1 — xx) on the left-hand side implies that the next-to- final entries in the 
symbol of /(x, x) contain the letter (1 — xx). 

In ref. [35], slightly simpler, but very similar, equations were analysed for a class of 
generalised ladder integrals. The analysis of ref. [35] can be adapted to the case of the 
four- loop integral J*^^^ and, as in ref. [35], the solution to the equation fl7.17p is uniquely 
determined by imposing single-valued behaviour on /. 
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First of all we note that any expression of the form h{x) — h{x) obeys the homogeneous 
equation and antisymmetry under the exchange of x and x and hence can be added to any 
solution of eq. fl7.17p . However, the conditions of single- valuedness, 

[disca; - discs] /(x, x) = , [disci„^ - disci_5;]/(a;, = , (7-18) 

and that and 1 are the only singular points, fix this ambiguity. 

Let us see how the ambiguity is fixed. Imagine that we have a single-valued solution 
and we try to add h{x) — h{x) to it so that it remains a single-valued solution. Then the 
conditions ( I7.18P on the discontinuites tell us that h can have no branch cuts at x = or 
X = 1. Since these are the only places that the integral has any singularities, we conclude 
it has no branch cuts at all. Since the only singularities of the integral are logarithmic 
branch points, h has no singularities at all and the only allowed possibility is that h is 
constant, which drops out of the combination h{x) — h{x). Thus there is indeed a unique 
single- valued solution to eq. fl7.17p . The argument we have just outlined is identical to the 
one used in ref. [31] to solve for the generalised ladders. 

A very direct way of obtaining the symbol of the single-valued solution to eq. (17.171) is 
to make an ansatz of weight eight from the five letters 

{x, 1 — X, X, 1 — X, 1 — xx} , (7-19) 

and impose integrability and the initial entry condition. Then imposing that the differential 
equation is satisfied directly at symbol level leads to a unique answer. 



7.3 An integral solution 



Now let us look at the differential equation (17.171) in detail and construct the single-valued 
solution. It will be convenient to organise the right-hand side of the differential equation 
(I7.17P according to symmetry under x 4-7- 1/x. We define 



E+(x,x) 



2 

E_(x,x) = 



X 



E(l -x,l 

E{l-x,l-x)- E[l 

Then the differential equation reads 

(1 — xx)xxdxdxf{x,x) = —(1 — xx)£'_(x,x) 
We may now split the equation (I7.2ip into two parts 



1 








1 - 


-)] 


X 




X/ - 


1 


1 - 




X 




-)] 

X/- 



'1 -|- xx)£'+(x, x) 



xxdxdsfa{,x,x) = -i?_(x,x) , 
'1 — xx)xxdxdxfb(yX, x) = — (1 -|- xx)E^{x, x) . 



(7.20) 



(7.21) 



(7.22) 
(7.23) 



Note that we may take both and ff, to be antisymmetrc under x O 1/x. 
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The equation f l7.22p is of exactly the same form as the equations considered in ref. 
Following the prescription given in ref. [33], Section 6.1, it is a simple matter to find a 
single-valued solution to the equation f l7.22p in terms of single-valued polylogs. We find 

fa{x, x) =L3^4_o — 2L4^-s,0 + -^5,2,0 + 2-^3^2,2,0 ^ 2L4_i^2,0 ~ -^4,2,0,0 — 2-^4^2,1,0 + -^5,0,0,0 

+ 2-/^5,i,o,o + 2-^v5^i 1 + 2L4 0,0 — 4:(^3{L^ — 2-^3 2 + 2-^4 + 3L4^i) (7.24) 

We now treat the equation fl7.23p for f^. Let us split it into two parts so that x) = 
fi{x,x) + f2{x,x), 

(1 - xx)xxd^dsfi{x,x) = -E+{x,x), 

(1 - xx)d^dsf2{x,x) = -E+{x,x). (7.25) 

We may write integral solutions 
and 

J2{x.x) = -Ml/x,l/x) = -l\tl\i^^^ (7.27) 

which obey the equations fl7.25p . 

It follows that the full function / is given by 

f{x, x) = fa{x, x) + /i(x, x) + /2(x, x) + h{x) - h{x) (7.28) 

for some holomorphic function h. We note that /(x, 1) — /a(x, 1) = h{x) — h{l). 

Now we examine the function /2 in more detail. Writing E+(t,i) = J2i^Wi(t)Hi,;'_{i) 
we find 

dt 

-jH^Xt)h,,{t,x) (7.29) 
where, for a word w made of the letters and 1, 



Iw {X-i 2^) 



^H^{^ = (-l)'^(G(i, w; x) - G(i, w- 1)) . (7.30) 



h t 

We may now calculate the symbol of /2. We note that 

df2{x,x) = dlogx'^ Hu,Xx)Iw[{x, x) - {x ^ x) . (7.31) 

i 

The symbol of is obtained recursively using 

S{Ux,x)) = SiHUx)) ® - {-ir'S{I^ix,x)) ® , (7.32) 

1 — xao 1 — xaQ 
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where w = QqW 



'. When w is the empty word /( logarithm, 



I{x, x) = log 



1 — XX 



(7.33) 



1 — X 



Using the relations f l7.31f7.32f7.33p we obtain the symbol of f2{x,x). One finds that the 
result does not obey the initial entry condition (i.e. the first letters in the symbol are 
not only of the form u = xx ot v = [1 — x){l — x)). However, the inital entry condition 
can be uniquely restored by adding the symbols of single-variable functions in the form 
<S{h2{x)) — S{h2{x)). Inverting x 1/x we may similarly treat fi{x,x) = — /2(l/x, l/x). 
Combining everything we obtain the symbol 



5(/(x,x)) =5(/,(x,x)+/2(x,x)+/i2(x)-/i2(x)-/2(l/x,l/x)-/i2(l/x)+/i2(l/x)). (7.34) 



The symbol obtained this way agrees with that obtained by imposing the differential equa- 
tion on an ansatz as described around eq. (17.191) . 

Given that single- valuedness uniquely determines the solution of the differential equa- 
tion f l7.17p one might suspect that we can use this property to give an explicit representation 
of the function h2{x). Indeed this is the case. The integral formula fl7.29p can, in principle, 
have discontinuities around any of the five divisors obtained by setting a letter from the 
set (I7.19P to zero. 

Let us consider the discontinuity of f2ix,x) 



The above expression vanishes due to the symmetry of Ej^ under x l/x and the anti- 
symmetry under x -H- x. The absence of such discontinuities is the reason that we split the 
original equation into two pieces, one for fa and one for fi,. 
Now let us consider the discontinuity around x = 1. We find 



The first term above again vanishes due to the symmetries of The second term will 
cancel against the corresponding term involving disci_fi?+(t, f) in the integrand when we 
take the combination [disci-^^ — disci_s]/2(x, x). The discontinuities at x = 1 and x = 1 of 
/2 therefore satisfy the single-valuedness conditions fl7.18p since does. 




(7.35) 




(7.36) 
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For the discontinuities at x = we find 



disc^/2(x,x) 



dt disc+ 



dt- 



tt 



dt 



dt 



1-tt 



discjE+(t,t) . (7.37) 



Now writing the t integral above as = — and using [disct — disct] E+{t,i) = we 
find 



[disc^ - discs] /2(x,x) 



dt 



dt 



tt 



discti?+(t, i) 



disc^ 



dt 



dt 



1 - tt 



E+it,i) 



+ (x O x) 

+ (x -H- x) 



(7.38) 



Thus /2(x,x) is not single- valued by itself since the above combination of discontinuities 
f l7.38p does not vanish. Note however that eq. f l7.38p is of the form A;(x) + A;(x), as necessary 
in order for it to be cancelled by adding a term of the form h2{x) — /i2(x) to f2{x,x). We 
now construct such a function h2{x). 
Let 

dt 



dt 



1-tt 



E+{t,t). 



Writing E+{t,t) = EiHn.At)H^'St) find 

dt 



hl{x) = - 



Y,H{wi-t) 



dt 



i-^/t 



H{wl;i). 



We can write 



A,^:(o,i);i) 



(7.39) 



(7.40) 



(7.41) 



where we have made explicit that w'^ is a word in the letters and 1 and d is the number 
of 1 letters. One can always rewrite this in terms of HPLs at argument t. Indeed we can 
recursively apply the formula 



/"* dr /"* dr 

G{la2,as...,a„-l)= / -^(f , as, . . . , a„; 1) - / — G(i, as, . . . , a„; 1) . (7.42) 

Jo ^ ~ Jo ^ 

to achieve this. Note that a^ G {0, 1} in the above formula. We also need 

G(l,0,;l) = (-l)^+iiJ,+i(t). (7.43) 

Once this has been done, one can use standard HPL relations to calculate the products 

H^^it)Gil/t,w',;l) (7.44) 

and perform the remaining integral from to x w.r.t. t. We thus obtain a function /ig 
whose discontinuity at x = is minus that of the x-dependent contribution to [disc^^ — 
discs] /2(a:,x). 
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In ref. [53] an explicit projection operator T was introduced which removes the dis- 
continuity at X = of a hnear combination of HPLs while preserving the discontinuity 
at a; = 1. The orthogonal projector (1 — J-") removes the discontinuity at x = 1 while 
preserving that at x = 0. We define 

h^ix) = (1 - 7)h\[x) . (7.45) 

Explicitly we find 

h2{x) = — — C6-f^2 + 17(3-^2 - -(2(3-^^2,0 ^C5-f^2,0 - (2(3-^^2,1 + -rC4-f^2,2 + 2C3-f/2,3 

21 19 17 3 

— C2-f^2,4 + — C4-f^2,0,0 + ^C4-f^2,l,0 + ^C4-f^2,l,l + 5C3-f^2,l,2 — C2-f^2,l,3 ~ 2^3-^^2,2,0 

1 3 

+ C3-f^2,2,l — C2-ff2,2,2 + -^C2H2,'i,Q — C2-f^2,3,l " 2*^3-^^2,0,0,0 " 3^3-^/2,1,0,0 — 3C3-f^2, 1,1,0 

1 11 

+ 4C3-f^2,l,l,l + C2-f^2, 1,2,0 ~ 2(^2 -f/2, 1,2,1 + --^2,1,4,0 + C2-f^2,2,l,0 + 2-^2,3,2,0 ~ 2-^^,4,0,0 

~ -f^2,4,l,0 + 2*52-^^2,1,0,0,0 ,1,1,0,0 + 2C2//2 ,1,1,1,0 ,1,1,3,0 ,1,2,2,0 ~ H2 ,1,3,0,0 

+ 2/^2,1,1,1,2,0 ~ -^2,1,3,1,0 + -f^2,2, 1,2,0 " 2-^^,2,2,0,0 ~ -^^2,2,2,1,0 + 2-^2,3,0,0,0 ~ -^^2, 3, 1,1,0 
~ -^2,1,1,2,0,0 + 2-^2,1,2,0,0,0 -f^2, 1,2, 1,0,0 ~ 2/^2,1,2,1,1,0 + -^2,2,1,0,0,0 + -^^2,2, 1,1,0,0 

+ ^2,1,1,1,0,0,0 . (7.46) 

Here the H functions are all implicitly evaluated at argument x. 

The contribution from fi{x,x) = — /2(l/x, 1/x) is made single- valued by inversion on 
X. So we define 

h{x) = h2ix)-h2il/x). (7.47) 
Finally we deduce that the combination 

/(x, x) = faix, x) + /i(x, x) + /2(x, x) + h{x) - h{x) (7.48) 

is single-valued and obeys the differential equation f l7.17p and hence describes the four- 
loop integral /(^) defined in equation ([Ip. The equations ([729]), fTTiH]) . fTriTj) and flTi8|l 
therefore explicitly define the function /. 

7.4 Expression in terms of multiple polylogarithms 

Now let us rewrite the integral form (17.291) . fl7.30p for /2(x,x) in terms of multiple poly- 
logarithms. We use the following generalisation of relation f l7.42p . 

G(i, 02, ... , an, z) ={G{\] y) - G{^] y))G{a2, . . . , a„; 2;) 

+ £(^-f)G(i,a3,...,a.;z) (7.49) 
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to recursively rewrite the appearing in the J^^(t, x) in eq. fl7.29p so that the t 

appears as the final argument. Note that in eq. fl7.49l) . the two terms involving an explicit 
appearance of l/a2 vanish in the case 02 = 0. The recursion begins with 

GC-;z) = hg{l-yz) = G{l;y). (7.50) 

The recursion allows us to write the products Hy^^(t)I^i,{t, x) as a sum of multiple polylog- 
arithms of the form G{w; t) where the weight vectors depend on x. Then we can perform 
the final integration dt/t to obtain an expression for /2 in terms of multiple polylogs. 
We may relate fi{x,x) directly to f2{x,x) since 

= hix^x) - 5^[ifo».(x) - Ho^XmHoy^'S^) - Ho^'Sl)] . (7.51) 

i 

For a practical scheme we express as a sum over H^-{t)Hyj'\t) and do the i integration. 
In any single term of the integrand of /2, the recursion (17.491) will lead to multiple poly- 
logarithms of the type G{. . . ,l/x;t). Next, we take the shuffle product with the second 
polylogarithm and integrate over t. 

In this raw form our result is not manifestly antisymmetric under x ■v^ x. Remarkably, 
in the sum over all terms only G(0, 1/x, . . .;x) remain. Upon rewriting 



r dt 

G (0, |,ai, . . . ,a„; = G(0; a;)G (|, ai, . . . , a„; a;) - / ^^(0; t)G'(ai 



, an, t) 



(7.52) 

we can use (I7.49P to swap G{l/x, . . . ; x) for (a sum over) G{. . . , x). Replacing the 
original two- variable polylogarithms by 1/2 themselves and 1/2 the x, x swapped version, 
we can obtain a manifestly antisymmetric form. The shuffle algebra is needed to remove 
zeroes from the rightmost position of the weight vectors and to bring the letters 1/x, 1/x 
to the left of all entries 1. Finally we rescale to argument 1. 

In analogy to the notation introduced for the Hard integral let us write 

G'3,2,1 = G(0,0,^,0,1,1;1) (7.53) 

etc. Collecting terms, we find 
l[%,{x,x) = (7.54) 

— -^2,2,4 + 2L2,3,3 — i^2,4,2 — 2^2,1,1,4 + 2L2, 1,2,3 — 21/2,1,3,2 + 2^2,1,4,1 + 2^2,2,1,3 — 2-^2,2,2,2 
~ 2X2,2,3,1 + 2^2,3,1,2 + 2L2,3,2,0 + 2Iv2,3,2,l — 2L2,4,l,o — 2L2,4,l,l — 2L3, 1,3,0 + 2L3,3,i,o 
~ 41/2,1,1,2,2 + 4-^2, 1,2, 2,1 + 4-^2,2,1,1,2 ~ 4L2,2,2,1,0 ~ 4-^2,2,2,1,1 ~ 4-^3,1^, 2,0 + 2-^3,2,1,0,0 
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+ 4^3,2,1,1,0 + Lq (-^/'i,2,4 + "^^1,3,3 ~ -^1,4,2 ~ 2^^1,1,1,4 + 2^^1,1,2,3 ~ 2/7^ ;^ 3 3 + 2H^^^^^^ 

+ 2i?]^ 2,1,3 " 2iyi 2,2,2 ~ 2i?i 2,3,1 + 2Hi ^ i 2 + 2^^1,3,2,1 ~ 2if^ 4 — 4:H^^^ -^^ 2,2 + 4^^i, 1,2,2,1 

+ 4i^l"2,l,l,2 - 4i^i~2,2,l,l) + 4i^l"2,5 - 4i^l"3,4 " 4i?i"4,3 + 4^^l"5,2 + 8-^1^1,1,5 " ^^1,1,2,4. 
+ 4^^1~1,4,2 - 8ifi"i 5 1 - 4:H[2,1A + 8i^l~2,3,2 + 4^^1~2,4,1 " 8i?i"3,i,3 " 4i/'i"4,i,2 - '^H{i 2,l 
+ 8-f^l~5,l,l + 8^^1^1,1,2,3 + 8^{^i, 1,3,2 - 8-f^l~l,2,3,l ~ ^^^,1,3,2,1 ~ 8-^1^2,1,1,3 + 8^1^2,3,1,1 
~ 8i/i 31 1 2 + 8i/i,3,2,i,i + C3 (8/^2,3 ~ 12/^3,2 — 121/2,1,2 + 12/^2,2,1 " 12/^3,1,0 — 16^3,1,1 

- 12Loi/i-i,2 + 12Loi/i-2,i + 16iJi;i_3 - I6ifi;3^i - i6i/i:i^i_2 + i6ifi;2,i,i) + 2L3C5 
+ Cl (-24L2 - 721/2" - 48iJi;i) + 

(—^4,2 — ^4,0,0 — 2L4,i,o — 2^4,1,1 — 4Z2,lC3 + 4LoC3^^2" + ^Li(3H2 

+ 3L2i?4 + 2Li^iH^ — ALiH^ + 4(^3if]^ 2 ^ 2Lo,o^i,3 + ^LqH^^ — 2LiHi^ — 12if^ 5 

+ 12C3-f^2^1 - -^0,0-^^2^2 + 2^0-f^2^3 ~ '^^1^2,3 ~ 6^2^4 + -^0-f^3^2 ~ '^LiH^2 ~ 4^0-f^4^1 

- 2LiH^i + 27/4^2 + ISi/g^i + 16C3^i"i,i — 2Lo,o^^i~i,2 + 6ivoi/i~i,3 — 12iJj^^ 4 — 2Lo,o^^i~2,i 

+ QLqH^2,2 ~ 8^1^2,3 + 2^0-f^l^3,l ~ 8ifi~3 2 + 4i/{^4,i + 2Lo,o^2^i,i + 2^0/^2^1,2 ~ 4/^2^^ 3 

- 2LqH2^2,i ~ 4^^2^2,2 + 4i/2^3,i - 6Loi/3;i,i + 4if3;i 2 + 12/^3^2,1 + 12i/4^i 1 

+ 4Loi?i"i,i,2 - 8i^i"i,i,3 - Si^i'i 2,2 - 4Loi?i~2,l,l + 8i^'i;2,2,l + 8i^i"3,i,i - ^2^^4^) + 

Gf (2X3,2 + 4X3,1,0 + 4X3,1,1 - 8(3^^2" - 6X2^3" - 2X0,0^/3" - 4Xi,ii?3- + 8LoH^ 
+ QLiH^ - 20H^ - 8C3H^^ + 8Loi/{;3 + AL^H^^ - 20H^^ + QLqH2,2 + 4Lii72"2 

- 1^2,3 + AL^Hi^ + 4LiiJ3-;i - IQHI2 - ^Hii + 4Loi/i-;i,2 - iQH^,i,^ + 4Lo^ri-;2,i 

- 16ifi"2,2 - 8i^'i^3,i - 4Loi^2^i,i - 8H2^^ 2 + 8^^3^i,i - 8i^i~i i 2 + 8ifi~2,i,i + 2L2H^) + 

Gt^ (2X3,2 - X4,0 - 2X4,1 + 2X3,1,0 - I2L2C3 + 16C3^^2~ - 4X2^^3" 

- 2U,oH^ + 8LqH^ + 8LiH^ - 20H^ + lQC3H^,i + 4:LqH^;^ - 16H^^^ + ALoH2^2 

- 12if2~3 - 12if3~2 + 4Loi/{^i 2 - 8i/i^i 3 - 8i/]^2,2 ~ 4Lo//2^1,l + 8i^2~2,l + 8^3^1,1) + 

(jr|^ (3LqH^ — 12H^ + 3LqHi 2 ~ 12i/i 3 + 6/^0^2 1 — ^'^^2 2 — ^'^^3 1 "I" ^LqH^ 

- 12iJ{;i 2 - 12i/i-;2,i) + Gf^ (2L3,o + 4L3,i - 8(3^^1" + 2/^0^3" " 4^ii^3" " 

+ ALoH^^2 - 8H^,3 + ^L^H2^, - 8H2^2 - 8H{^i,2 + 8^^2:1,1) + Gf^ + ^^3,0 + 4L3,i 

- leCai^r - ^L^H^ + ^H^,3 + ^H2^2 + ^Hl, - 4Lo//i;i,i + 8i7i-2,i + 8H2^,^^) + 

(-2L4 + 2L3,o + 16C3^r + 41^0^^1:2 - ^Hl, - 8H2^2) + 
Gf {-AH^ - AH^2 - 8^2:1 - 8i^i;i,i) + Gt^ (-6X2,1 + ?>LqH2 + QUH2 - 12H^ 

- 12H^2 - 12i^2Ti) + GI2 (~2^3 - 8X2,1 + 2L0H2 + 8L1H2 - AH^ - 8H^2 

- 8H2^, + 8/Ji-i,i) + ^3+^^^ (4X3 + AL0H2 - IQH^ - 8H^^2) + (-4X3 - 8X2,1 
+ 2L0H2 + 8L^H2 - AH^2 - ^H2^i + 8i/i:i,i) + Gt^^^ (4X3 - 8H^ - 8H^^2) + 

2 (4^3 + 4X2,1 - ^UH2 - AH^ + 8H2^,) + G+^_^_^ (4Loi/2" - 8H^) + 
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4Gi 



{2,2 ~ ^'^tsA ~ ^'^tl,2 ~ 8^^2,1 ~ 6*^41,1 



'2,1,1,2 --2,2,1,1^ - ^ - I2GJ3 - 6Gt^ - AGl + 12^^ 3 

(-16^4 - - 12G^, - 8GI + 8GJ^ 3 + SGJ^^^ + 8GlJ H^,+ 

+ + I2GJ3 + 4GJ^ - 16Gl^^ - 12Gj^ 3 - 12Gj3^^ - 16Gj^ 3 

8G^„ , - 12G 



- 16G 
-8G 



3,2,2 ---3,3,1 1,2 ~ ^2G~;2,1 ~ 1,1 + 1,1,3 + 1,2,2 " ^^^2,2,1 



2,3,1,1 



!l,l,2 - 8^^2,1,1) ^1 



^ (2/73^2 — -^4,0 — 2/74^1 + 1^3^0,0 + 2L3^i_o — I2L2C3 + — 4:L2H^ 

+ 8LiiJ+ - leCsi/iti - ^H+, + 4LoiJ3+ - 4if3+ - 1QH+, - 4LoiJ+ ,2 + 8Hli,3 
+ 8^1^2,2 + 4^oi^2ti,i - 8Ht,2,i - 8Hli,i) + 

G7^ (2X3,0 + 4L3,i + 4L1C3 - 6LoH^ - ALiHt + 2QHt - ^LoH^^ + IS^ts 

- ALoHl^ + + Si^sti + 8i^i+i,2 - ^H^i^) + (^4 + ^3,0 + 4X3,1 + 8L1C3 

- 2L^Ht - AL^Ht + mt + 4/7+ + AH^^ + 4i:oi^iti,i - 8i/i+2,i " ^^1,1) + 
<^2~,i,i (-2I4 + 2X3,0 - 8L1C3 - ^LoHt + IQHt - + ^H^^s + 8^^212) + 
G[ 1 (-3^2,0 - 6L2,i + SLoHt + eLi//^") + (-2L3 - 2^2,0 - 8^2,1+ 
2LoH^ + 8LiH^ - 8iy+2 - 8i72ti - 8Ht^i,i) + G^^^^ (4L3 - 2^2,0 + 8H+2) + 

(—4-^3 — 2^2,0 — 8-^2,1 — i^o,o,o — 8C3 + 2L0H2 + 8L1H2 — ^2H^2 

- 8H+, - 8i/+ ,1) + Gr^_^ (4L3 + 2Lo,o,o + I6C3 + 8H+2) + G^r^ ^ (^^3 + ^^2,0 

+ 4L2,i + ^0,0,0 + 8C3 - 2L0H+ - AL,H+ + 8H+2) + G^;, , ^ (-2^2,0 - 2Lo,o,o - I6C3) + 
Gr^ {4H^ + 8ff+ ) + (3I2 + 12i^+i) + ^1,1 (-6Z2 + 12H^) + ^3^3 (8I2 + 2Zc 

- AH+ + 16/7+ ) + ^3-^^^ (-8I2 - 4Zo,o + 8H+) + G7^^^ {-8L2 - 2Zo,o + 8Ht 

- 8H+,) + ^3^^^^^ (4Lo,o + 8H+) + Gr^ (IOL2 + 4Lo,o - + ) + 
'^2\3,i (-4X2 - 4X0,0) + G7^^^^ (-8X2 - 4X0,0 + - 8H+,) + G2:;2,i,i ("4^2 + 8/^2+) + 

-10G7^ -8Gr^ 

3G 

3,^,2 3,3,i 

+ 4G 
f-lOG 



^0,0 



2,1,3 V --2 - 4X0,0 + - 8H+,) + G^:;^ 2,1 ^^^O'O + 81^2') + ^2:;,^^ 2 (^^2 + 4Xo,o) + 

2,5 -^3,4 --4,3 ' 10^2.1,4 + 8^*2^2,3 + 8^2,3,2 + 4^Ul + 8^11,3 



6G7 



3,3,1 ' -4,1,2 ■ -4,2,1 -8Gr^^^^3-8Gr „„-4Gr..^,-4G 



2,3,2 

^^2,1,2,2 



2,1,3,1 



2,2,1,2 



2,3,1,1 - 8^11,1,2 - 8^11,2,1 - , . + 4G 



4,1,1,1 



- 4Gr 

2,1,1,1,2 2,1,2,1,1' 



. , - lOGr^ - 6Gr,^ - 2Gr^ + 8G: 
2,5 3,4 4,3 5,2 



2,1,4 ' 6*^2,2,3 + ^^2,3,2 + 8<^3,l,3 



+ 8G 

+ 4G 



3,2,2 ' ^^3,3,1 ' 6^4,1,2 + 6^4,2,1 + 4*^5,1,1 ~ 4^2,1,1,3 ~ ^^2,1,2,2 + 4^^2,2,2,1 



4G: 



2,3,1,1 " 4Gg^^^^ 2 + 4<^3,2,l,l) ^1 
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V 


Analytic 


r LtLO i A 





n 1 
U.i 


n o 


iOD. ( oo 


iOD. ( oo 


A Oq 7 

4. ye- ( 


0.2 


0.3 


116.962 


116.962 


5.9e-8 


0.3 


0.1 


110.366 


110.366 


2.8e-7 


0.4 


0.5 


84.2632 


84.2632 


1.4e-7 


0.5 


0.6 


75.2575 


75.2575 


1.4e-7 


0.6 


0.2 


78.3720 


78.3720 


3.7e-8 


0.7 


0.3 


70.7417 


70.7417 


6.8e-8 


0.8 


0.9 


58.6362 


58.6363 


1.4e-7 


0.9 


0.5 


60.1295 


60.1295 


l.le-7 



Table 4: Numerical comparison of the analytic result for ^^3X24 /''^■'(xi, X2, X3, X4) against 
FIESTA for several values of the conformal cross ratios. 



20Gr6 + 20^3^, + 12Gr^ + ^G^^ - 20G- - 16Gr^^^ - 12^^^^ 3 - 12Gr^^^ 
- 20^3-,^, - 16^3^^ 3 - 16^3-3^^ - 8^3-^^^ - 12Gr_^ 3 - 12Gr^_^ - 12^73^^ - 4G^7^^^ 

~ ^^§,2,1 + ^6^2,1,1,4 + ^2^^1,2,3 + ^2G^^i,3,2 + 8^*^2,1,3 + ^^2^2,2,2 " ^^M,2,l ~ 

+ 16G3-,,3 + 16^3-^^,^, + 8^3-^ 3_^ + 8^3-^^^^^ - 8^3-3^^^^ + 12Gr^^^_^ + 12Gr^^^^^ + 8Gr^^^^^ 

~ ^^2;i,l,l,3 ~ ^^2,1,1,2,2 + ^^M,2,2,l + ^^^^l.S,!,! " ^'^UX1,2 + ^^3!l,2,l,l 

7.5 Numerical consistency tests for I^^^ 

In order to check the correctness of the result from the previous section, we evaluated 
I^^^ numerically and compared it to a direct numerical evaluation of the coordinate space 
integral using FIESTA. In detail, we evaluate the conformally-invariant function f{u,v) = 
^13^24 -^^^H^i' ^2; 2:3, X4) by first applying a conformal transformation to send 0:4 to infinity, 
the integral takes the simplified form, 

"i" •^13-^24 ^14;23 - ^ 2 ^2 ^2 ^2 ^2 ^2 ^2 ^2 ^2 ' ^i-OO) 



and then using the remaining freedom to fix x\.^ = 1 so that and f = X23. In com- 

parison with the two 3-loop integrals, the extra loop in this case yields a moderately more 
cumbersome numerical evaluation. As such, we modify the setup for the 3-loop examples 
slightly and only perform 5 x 10^ integral evaluations. We nevertheless obtain about 5 dig- 
its of precision, and excellent agreement with the analytic function at 40 different points. 
See Table S] for an illustrative sample of points. 
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8 Conclusions 



Recent years have seen a lot of advances in the analytic computation of Feynman integrals 
contributing to the perturbative expansion of physical observables. In particular, a more 
solid understanding of the mathematics underlying the leading singularities and the classes 
of functions that appear at low loop orders have opened up new ways of evaluating multi- 
scale multi-loop Feynman integrals analytically. 

In this paper we applied some of these new mathematical techniques to the computation 
of the two so far unknown integrals appearing in the three-loop four-point stress-tensor 
correlator in A/" = 4 SYM, and even a first integral occurring in the planar four-loop 
contribution to the same function. The computation was made possible by postulating 
that these integrals can be written as a sum over all the leading singularities (defined 
as the residues at the global poles of the loop integrand), each leading singularity being 
multiplied by a pure transcendental function that can be written as a Q-linear combination 
of single-valued multiple polylogarithms in one complex variable. After a suitable choice 
was made for the entries that can appear in the symbols of these functions, the coefficients 
can easily be fixed by matching to some asymptotic expansions of the integrals in the limit 
where one of the cross ratios vanishes. In all cases we were able to integrate the symbols 
obtained form this procedure to a unique polylogarithmic function, thus completing the 
analytic computation of the the three-loop four-point stress-tensor correlator in A/" = 4 
SYM. While for the Easy integral the space of polylogarithmic function is completely 
classified in the mathematical literature, new classes of multiple polylogarithms appear in 
the analytic results for the Hard integral and the four-loop integral we considered. 

One might wonder, given that the Hard integral function H^'^^ involves genuine two- 
variable functions, whether there could have been a similar contribution to the Easy inte- 
gral, compatible with all asymptotic limits. Indeed there does exist a symbol of a single- 
valued function, not expressible in terms of SVHPLs alone, which evades all constraints 
from the asymptotic limits. In other words the function is power suppressed in all limits, 
possibly up to terms proportional to zeta values. However, the evidence we have presented 
(in particular the numerical checks) strongly suggests that such a contribution is absent 
and therefore the Easy integral is expressible in terms of SVHPLs only. 

We emphasise that the techniques we used for the computation are not limited to 
the rather special setting of the = 4 model. First, by sending a point to infinity a 
conformal four-point integral becomes a near generic three-point integral. Such integrals 
appear as master integrals for phenomenologically relevant processes, like for example the 
quantum corrections to the decay of a heavy particle into two massive particles. Second, 
the conformal integrals we calculated have the structure ^ Ri Fi (so residue times pure 
function) that is also observed for integrals contributing to on-shell amplitudes. However, 
we believe that this is in fact a common feature of large classes of Feynman integrals (if 
not all) and one purpose of this work is to advocate our combination of techniques as a 
means of solving many other diagrams. 

Further increasing the loop-order or the number of points might eventually hamper our 
prospects of success. Indeed, beyond problems of merely combinatorial nature there are also 
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more fundamental issues, for example to what extent multiple polylogarithms exhaust the 
function spaces. It is anticipated in ref. pU|IM] that elliptic integrals will eventually appear 
in higher-point on-shell amplitudes. Via the correlator /amphtude duality this observation 
will eventually carry over to our setting. Nevertheless, some papers [6T1|62] also hint at a 
more direct albeit related way of evaluating loop-integrals by casting them into a '(i log- 
form', which should have a counterpart for off-shell correlators. 
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A Asymptotic expansions of the Easy and Hard inte- 
grals 

In this appendix we collect the asymptotic expansions of the different orientations of the 
Easy and Hard integrals in terms of harmonic polylogarithms. The results for -E'i4;23 and 
Hi2;34 were already presented in Section [31 The results for the other orientations are given 
below. 



^13^24 -^12;34 



log U 
log^ U 

1 



1 

3x 



{A.l] 



X' 



^2if2,2 + -f^2,l,l + 2i/i,3 + Hii^ 



2x 
\ogu 



4i?2,2 — 3if2,i,i — 4-^1,3 — i?i,2,i — 4i/i 



1,2 



— ( — 16i/3,2 — 8//34,i — 16/^2,3 — 8i/2,l,2 " 8i/i,4 + 4 -ffl,3,l 



— 4i^i 2,2 — -f^l,2,l,l — 4iJi 1 3 + 2 i/l, 1,2,1 — -f^l,l,l,2 

+ - {SH^ 2 + 5-^311 + 8iJ2,3 + -f^2,2,l + 6-f/'2,l,2 + 4 i^i 4 — -f^l,3,l 

X V 

+ 5iJi_2,2 + -f^l,2,l,l + 4iJi,i,3 — 2 -f^l, 1,2,1 + -f^l,l,l,2 
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+ ^ (4C3i/l,2 + 2(3^1,1,1 + 32/^4,2 + 16i/4,l,l + 32 i/3,3 + 16i73,l,2 

+ 167^2,4 — 8iJ2,3,l + 87^2,2,2 + 2 i/2,2,1,1 + 8iJ2,l,3 ~ 47^2, 1,2,1 ,1,1,2 

— 8ifi,4,i + 4 -ffl,3,2 + 4-^1,3,1,1 + 4i7i,2,3 — 2-f/'i,2,2,l + 2-^1,2,1,2 — 4 -ffl, 1,3,1 

+ -f^l,l,2,l,l — -^^1,1,1,2,1 j + ~ ( ~ 4^3-^^2,1 — 6C3-f?^l,2 — 2^3-^^^1,1,1 — 16 -^^4,2 

— 10-ff4,i,i — 16-ff3,3 — 10-ff3,i,2 — 8-ff2,4 + 4-^2,3,1 ~ 8 -f^2,2,2 " 2-^^2,2,1,1 

— 6-f/'2,l,3 + 4-f^2, 1,2,1 — 2-f^2,l,l,2 + 2 -??1,4,1 " 6-f/'i,3,2 — 4-ffi,3,i,i — 6-?fi,2,3 
2-f^l,2,2,l ~ 2-^1,2,1,2 + 4-^^1,1,3,1 — -^1,1,2,1,1 + -^1, 1,1, 2,1 ~ 8(^3 -f^3 



+ 20C5^i) + 0{u) 



^13^24 -E'l3;24 



logM 



X 



{H2,2,1 — -f^2,l,2 + -^1,3,1 ~ -^1,2,1,1 — -^1,1,3 + -f^l, 1,2,1 " 6 C^H-} 



+ - (4(3-^^2,1 — 2(3-^^1,2 — 2-^/3 2.1 + 2H312 — 2-/^2,3,l + -f^2,2,l,l + 2-^2,1,3 

— 2-ff2, 1,2,1 + -f^2, 1,1,2 — 4-f/^l,4,l + 3-^1,3,1,1 + -f/^1,2,1,2 + 4-f^l,l,4 " 2-^11,1,3,1 

- i^l,l,2,2 - ^^1,1,2,1,1 - i^l,l,l,3 + ^1,1,1,2,1 + 12(3^^3) + ^(m) , (A.2) 

xt3X^4 ^13;24 = V « (2i/2,l " ^1,2 " i^^l,l,l) + ^(^^^ (^2,1 " ^^3) ] (A.3) 



+ log^ u 

1 

+ 



— ^ — 4-ff3,i — 2-^2,2 + 2-^1,3 + 2-f/'i,2,l 



1,1,2 



^ ~ 2-ff3 1 — H2,2 — -f^2,l,l — -f^l,3 + -f^l,2,l + 4-^4 ) 

(1 — XjX V / J 

+ log U [-^ (l6i/3,2 + 8i/3,l,l + 8iJ2,3 - 8i/2,2,l " 4ifi,4 - 12ifi,3,l 

(41/4,1+41/3,2 



— 4-f/'i 2,2 + 2-f/i 2,1,1 — 4-f/'i 1 3 — 2-f/'i 1 1 2) + 7:; r— 

/ (1 — XjX 

+ 6-f/3,i,i + 4-f/2,3 + 2-f/2,l,2 + 8-f/i,4 — 4-f/i,3,i — 2-f/i,2,2 — 2-f/i, 2,1,1 



^ (32(31/2,1 - 16(3^^1,2 - 16C3i^i,i,i 
... , .... .... _24i/2,4 



— 2-f/i,i,3 + 2-f/i,i,2,l — 20-f/5y J , ^2 

+ 64i/5,l - 32i/4,2 - 32i/4,l,l - 24i/3,3 + 16i/3,2,l - 16//3,1,2 - 

+ 401/2,3,1 - 41/2,2,1,1 - 8//2,l,3 + 4/^2,1,1,2 + 40/^1,4,1 + 4//i,3,2 

— 8//i,3,i,i — 4//i,2,3 + 4//i,2,2,l — 4//i,2,l,2 + 8//i,i,i,3j 
+ ^ , f 16(3^^2,1 - 4//4,2 - 12//4,1,1 - 4//3,3 " 12//3,2,1 " 8//3,l,2 

(1 — x)x \ 

— 12//2,4 + 4//2,3,i + 2//2,2,l,l — 4//2, 1,2,1 + 2//2,i,i,2 — 20//i,5 + 4//l,4,l 
+ 4//i,3,2 + 6//i,3,i,i + 4//i,2,3 + 2//i,2,l,2 + 8//i,i,4 — 4//i, 1,3,1 — 2//i, 1,2,2 

— 2//i,i,2,i,i - 2//i,i,i,3 + 2//i,i,i,2,i - 16(3^/3 + 40//6) + , 
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2^13^24 Hii.23 = — f — ('^H2,l — Hi 2 + 2ifi i i) — 27^2,1 " -f^l.2 " 2i^l,l,l — 



+ log^ u 

X \ 



— (^^3 1 + 7^2,2 — -f^l,3 + 2i/i,2,l 

-f^l,3 + -f^l,2,l + Hi 1^2 + -f^4 j 



-f^2.2 



2-f^l,l,2 



(A.4) 



r4 / 

+ log n — 4if3,2 - 4i73.1,l + 2i/2,3 + 4iJ2,l,2 - Hi^4 + 2ifl,3,l + 4iJi,2,2 
- 2iJi,2,l,l + m,l,3 + 2/^1,1,1,2) + ^ (2i/4,l + 4i/3,2 - 2if3,l,l + 4i/2,3 

+ 2iJ2,i,2 + 5-^1,4 + 2_?7i,3,i + 4i7i,2,2 — 2_ffi,2,i,i + 4ifi,i,3 + 2_ffi,i,i,2 



+ 5i/5 



(4(3^^2,1 - 2C3^i,2 + 4C3i^i,i,i + 8H,^i - 4/^4,2 + 8/^4,1,1 



a; 



— 6if3,3 + 4if3,2,l — 47/3,1,2 — 3if2,4 + 2i/2,3,l — 4if2,2,2 + 2i^2, 2,1,1 



— 6if2,l,3 — 2-^2,1,1,2 + 2i^i,4,i — 3-ffi,3,2 + 4ifi,3,i,i — 5/^1,2,3 + 2-f/'i, 2,2,1 

\ 8 / 

1,2,1,2 — 4i/i,i,4 — 2ifi, 1,2,2 — 2i/i,i,i,3 j ~^ ~y 4^3-^^2,1 — 2^3-^^1,2 

-f^l 11 + 37/4 2 — 2iJ4,i,i + 37/3,3 — 2iJ3,2,l + 4iJ2,4 + 2iJ2,2,2 



- 2i/ 

- 4C3i/l,l,l 
+ 2/^2,1,3 + 5//i,5 

+ 2/^1,1,2,2 + 2ifi, 1,1,3 



37/4,2 — ^-04 

3-ffi,3,2 — 2/71,3,1,1 

2(3/73 + ST/e) + 0(w 



+ 3-f/'l,2,3 — 2ifi,2,2,l + 4-^1,1,4 



B An integral formula for the Hard integral 

We want to find an integral formula for pure functions which involve x — x in the symbol 
as well We are interested in single-valued functions, i.e. ones obeying 

the constraints on the discontinuities, 

[disca; - disc5;]/(x, x) = , [disci-^^. - disci_5;]/(x, x) = . (B.l) 

and with no other discontinuities. 

It will be sufficient for us to consider functions whose symbols have final letters drawn 
from a restricted set of letters, 



S{F) = S{X) ® - + S{Y) ® + S{Z) ®{x-x). 
X 1 — X 



(B.2) 



where X, Y, Z are single- valued functions of x, x. 

We will suppose also that the function F obeys -F(x, x) = 0, as required to remove the 
poles at X = X present in the leading singularities of the conformal integrals. We therefore 
take Z(x, x) = also. If F has a definite parity under x O x then X and Y have the 
opposite parity while Z has the same parity. 

The functions X, Y and Z are not independent of each other. Integrability (i.e. d'^F = 
0) imposes the following restrictions, 

A rflog§ + A dlog^ § + dZArflog(x-x) = 0. (B.3) 



X 



1 — X 
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We may then define the derivative of F w.r.t x to be 



dr,F{x, x) = — 

X 



Y Z 
- + 



so that 



F{x, x) 



dt 



1 — X X ~ X 

X{t,x) Y{t,x) Z{t,x) 



(B.4) 
(B.5) 

(B.6) 
(B.7) 



t 1-t t-x 

A trivial example is the Bloch-Wigner dilogarithm function, defined via, 

F2ix,x) = log XX (ifi(x) - Hi{x)) - 2{H2{x) - H^ix)) . 
It has a symbol of the form flB.2p where 

Xi = log(l — x)(l — x), Yi = — log XX Zi = . 
Thus we can write the integral formula (IB. 51) for F2. 

B.l Limits 

We want to be able to calculate the limits of the functions to compare with the asymptotic 
expressions obtained in Section [31 The formula (]B.5|) allows us to calculate the limit x — 
(which means dropping any power suppressed terms in this limit). We may commute the 
limit and integration 



lim Fix, x) 



dt lim 



X{t,x) Y{t,x) Z{t,x) 
t 1-t ^ 



t 



(B.8) 



In the second and third terms one may also set the lower limit of integration to zero, 
directly. In the first one should take care that contributions from X{t,x) which do not 
vanish as t — ?■ produce extra logarithms of x, beyond those explicitly appearing in the 
limit of X, as the lower limit approaches zero. 



B.2 First non-trivial example (weight three) 

The first example of a single-valued function whose symbol involves x — x is at weight 
three [30]. There is exactly one such function at this weight, i.e. all single- valued functions 
can be written in terms of this one and single-valued functions constructed from single- 
variable HPLs with arguments x and x only. It obeys F3(x, x) = — F3(x,x). The symbol 
takes the form f lB.2p with 



X2 = -log(xx)(/7i(x) + i7i(x)) + ^(iJi(x) + i7l(x))^ 
Y2 = ~ log^(xx) + log(xx)(iJi(x) + Hi{x)) , 

Z2 = 2 log XX (i7i(x) - H^{x)) - A{H2{x) - H2{x)) . (B.9) 
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Note that X2, Y2 and Z2 are single- valued and that Z2 is proportional to the Bloch-Wigner 
dilogarithm (it is the only antisymmetric weight-two single-valued function so it had to 
be). They obey the integrability condition flB.3l) so we can write the integral formula flB.5|) 
to define the function F3. 

We have constructed a single-valued function with a given symbol, but in fact this 
function is uniquely defined since there is no antisymmetric function of weight one which 
is single- valued which could be multiplied by C2 and added to our result. Moreover, since 
it is antisymmetric in x and x, we cannot add a constant term proportional to Cs- 

Looking at the limit x — )■ we find, following the discussion above, 

limF3(x,x) =^\o<^ xHx[x) ^\o<gx{H2[x) + HiJx) - Hii{x)) 

- 3H3{x) - iJi,2(x) + H2,o{x) + H2,i{x) + i/i,o,o(a;) - i/i,i,o(a;) (B.IO) 

Starting from the original symbol for F3 and taking the limit a; — )■ we see that the above 
formula indeed correctly captures the limit. 

B.3 Weight five example 

We now give an example directly analogous to the weight-three example above but at 
weight five. The example we are interested in is symmetric F5(x, x) = ^5(0;, x). It has a 
symbol of the canonical form (1B.2P with 

X4^{x,x) = (>Co,0,l,l — -^1, 1,0,0 ~ -^0, 1,1,1 + -^1, 1,1,0) ; 

Y4^{x,x) = (£0,0,0,1 — -^1,0,0,0 — -^0,0,1,1 + -^^1, 1,0,0) ; 

Zi{x,x) = (£0,0,1,1 + £1,1,0,0 — £0,1,1,0 — £1,0,0,1) • (B-11) 

The above functions are single- valued and obey the integrability condition and therefore 
define a single- valued function of two variables of weight five via the integral formula. 
Taking the limit x — > we find 

limF5(x,x) =Hi^iHofifi + (i^i,i,o - i/i,i,i)/fo,o 

x-j-O 

+ {—H-s^i + -f/^2,1,1 + -f^i, 1,0,0 — -f^i,i,i,o)-f^o 

— Hi^4 — i^2,3 + 2iJ4,i -|- -f^l,3,l — -^3,1,0 ~ -^3,1,1 
+ H2 

,1,1,0 ,1,0,0,0 " ,1,1,0,0 + 2f/i,iC3. (B.12) 

This formula correctly captures the limit taken directly on the symbol of F^. This weight- 
five function plays a role in the construction of the Hard integral. 

B.4 The function iif^"^ from the Hard integral 

The function H^""^ from the Hard integral is a weight-six symmetric function obeying the 
condition H^°^{x^x) = 0. The symbol of H^""^ is known but is not of the form (lB.2p . 
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However, we can use shuffle relations to rewrite the symbol in terms of logarithms of u and 
V and functions which end with our preferred set of letters. We find the symbol can be 
represented by a function of the form 

H^"'\l-x,l-x) 

= (2i7o,o(«) + ^Ho{u)H^{v) + 8ifo,o(^))(>Co,o,i,i + -^1,1,0,0 " -^0,1,1,0 - C^flfi,i) 
- SF^iHoiu) + 2Ho{v)) + Fe . (B.13) 

Here F5 is the weight-five function defined in Section IB.3[ The function Fq is now one 
whose symbol is of the form (]B.2p . where the functions X5, and Z5 take the form 

X5 = 20£o,o,o,i,i + 12£o,o,i,i,o ~ 32£o,o, 1,1,1 ~ 8-^0,1,0,1,1 ~ 12£o, 1,1,0,0 8-^0,1,1,0,1 
+ 16i2o,i,i,i,i — 8£i^o,o,i,i + 8£i^o,i,i,o — 20£ij,o,o,o + 8£ij,o,o,i + 8£i, 1^0,1,0 



+ 32/:i,i,i,o,o - 16/:i,i,i,i,o - 16/:i,iC3 , (B.14) 
Y5 = 20£o,o,o,o,i — 32£o,o,o,i,i ~ 8£o,o,i,i,o + 16>Co,o, 1,1,1 ~ 8£o,i,o,o,i + -^0,1,1,0,0 

— 20£i,o,o,o,o + 8£i^o,o,i,o + 16>Ci,o,o,i,i + 8£i,o,i,o,o + 32£i^i^o,o,o — 16>Ci, 1,0,0,1 

- 16£i,i,i,o,o - 16/:i,oC3 + 64/:i,iC3 • (B.15) 
= 32F5 . (B.16) 



Note that the ^3 terms have been chosen in such a way the the functions X5, Y5 and 
obey the integrability condition (IB.3p . The integral formula for Fg based on the above 
functions will give a single- valued function with the correct symbol, i.e. one such that if^"-* 
defined in eq. (IB.13P has the correct symbol and is single-valued. 
We recall that the Hard integral takes the form 

H^-){l-x,l-x) ^ 

{x — xY {1 — xx){x — x) 

Calculating the limit x — j- we find that H^"''' reproduces the terms proportional to l/x"^ 
in the limit exactly, including the zeta terms. Note that in this limit the contributions of 
and H^^'^ are distinguishable since the harmonic polylogarithms come with different 
powers of x. Since there are no functions of weight four or lower which are symmetric in x 
and X and which vanish at x = x and which vanish in the limit x — ?■ 0, we conclude that 
i/^*^) defined in eq. (1B.13P is indeed the function. Comparing numerically with the formula 
obtained in Section we indeed find agreement to at least five significant figures. 

C A symbol-level solution of the four-loop differential 
equation 

In this appendix we sketch an alternative approach to the evaluation of the four-loop inte- 
gral. Mor precisely, we will show how the function J*^^-' can be determined using symbols and 



H 



14;23 



X' 



13 



4 
X24 
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the coproduct on multiple polylogarithms. We start from the differential equation fl7.17p . 
which we recall here for convenience, 

d^dsf{x,x) = _ Ei{x,x) - E2{x,x) , (C.l) 

(1 — xxjxx (1 — xx) 

where we used the abbreviations Ei{x, x) = E{1 — x,l — x) and E2{x, x) = E{1 — 1 — 
1/x). We now act with the symbol map S on the differential equation, and we get 

d^d,S[hx,x)] = ^——S[E,{x,x)] - -J—S[E2{x,x)] , (C.2) 

(1 — XXjXX (1 — XX) 

where the differential operators act on tensors only in the last entry, e.g., 

d^[ai (g) ...» a J = [(9^ log a„] ai (g) . . . (g) a„_i , (C.3) 
and similarly for (9^. It is easy to see that the tensor 

= S[Ei{x,x)] (g) ( 1 - — ) (g) {xx) +S[E2{x,x)] (g) {1-xx) (g) {xx) (C.4) 

solves the equation flC.2|) . However, 5*1 is not integrable in the pair of entries (6,7), and 
so Si is not yet the symbol of a solution of the differential equation. In order to obtain 
an integrable solution, we need to add a solution to the homogeneous equation associated 
to eq. (lC.2p . The homogeneous solution can easily be obtained by writing down the most 
general tensor 5*2 with entries drawn from the set {x, x,l — x,l — x,l — xx} that has the 
correct symmetries and satisfies the first entry condition and 

d.d^S2 = . (C.5) 

In addition, we may assume that 5*2 satisfies the integrability condition in all factors of 
the tensor product except for the pair of entries (6, 7), because 5*1 satisfies this condition 
as well. The symbol of the solution of the differential equation is then given by Si + S2, 
subject to the constraint that the sum is integrable. It turns out that there is a unique 
solution, which can be written in the schematic form 



S[f{x, x)] = s7 (g) u (g) u + So (g) (g) w + s7 (g) (g)- + st(g)-(g)u 

^■'^ ' ^ ^ l-X X X 

+ st®U®- + st® (g)-U + st(g)f(g)- + Sg(g)-(g)- 

+ Sg (g) (1 - U) (g) U , 



(C.6) 



where sf are (integrable) tensor that have all their entries drawn from the set {x, x, 1 — 
x, 1 — x} and the superscript refers to the parity under an exchange of x and x. 

The form (1C.6P of the symbol of f{x,x) allows us to make the following more refined 
ansatz: as the sf are symbols of SVHPLs, and using the fact that the symbol is the 
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maximal iteration of the coproduct, we conclude that there are linear combinations f^{x, x) 
of SVHPLs of weight six (including products of zeta values and SVHPLs of lower weight) 
such that S[f^{x, x)] = sf and 



^6,1,1 [/(a;, x)] = fi (x, x) ®\ogU® log U + f2 {x, x) ®\ogV ® log U 

+ f^ix, x) ® log (g) log - + ftix, x) ® log - (g) logu 

\ — X X X 



+ ftix, x) ®\ogu® log - + faix, x) ® log \ — - ® \ogu (C-7) 

X 1 — X 

+ /^(x, x) ®\ogv ®i log - + /^"(x, x) (g) log - ® log - 

+ /9~(x, x) ® log(l —U)<®U. 

The coefficients of the terms proportional to zeta values and SVHPLs of lower weight 
(which were not captured by the symbol) can easy be fixed by appealing to the differential 
equation, written in the fornix 

(id®9^®aj.)A6,i,i[/(x,x)] = ^——Ei{x,x)®mi-—^—-E2{x,x)®mi. (C.8) 

(1 — xxjxx (1 — xx) 

The expression (1C.7P has the advantage that it captures more information about the func- 
tion /(x, x) than the symbol alone. In particular, we can use eq. (IC.7P to derive an iterated 
integral representation for /(x, x) with respect to x only. To see how this works, first note 
that there must be functions y4''=(x, x), that are respectively even and odd under an ex- 
change of X and x, such that 

^ X 
A7,i[/(x,x)] = y4"(x,x) (g) logu + y4+(x,x) ®\og- . (C.9) 

with 

X 

A6,i[A"(x, x)] = /f (x, x) ®\ogu + /2~(x, x) (g) logf + f^{x, x) ® log - 

1 — X 

+ ^) ® log + /9~(a;, x) ® log(l - u) , 



I — X 



(C.IO) 



A6,i[v4+(x,x)] = /a (x,x) ® log + fbi^^^) ® log« 



X 



+ (x, x) ® log V + /g (x, x) (g) log - . 



X 

The (6,1) component of the coproduct of A+(x,x) does not involve log(l — u), and so 
it can entirely be expressed in terms of SVHPLs. We can thus easily obtain the result for 
A~^{x, x) by writing down the most general linear combination of SVHPLs of weight seven 
that are even under an exchange of x and x and fix the coefficients by requiring the (6,1) 
component of the coproduct of the linear combination to agree with eq. (IC.lOp . In this 



''We stress that differential operators act in the last factor of the coproduct, just like for the symbol. 
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way we can fix A^{x, x) up to zeta values of weiglit seven (wliicli are integration constants 
of the original differential equation). 

The coproduct of A~{x,x), however, does involve log(l — u), and so it cannot be 
expressed in terms of SVHPLs alone. We can nevertheless derive a first-order differential 
equation for A~{x,x). We find 

9^v4"(x,x) = - [/r(x,x) + fl{x,x)] - [f2{x,x) +/6+(x,x)] 



1 — XX 

= K{x, x) . 
The solution to this equation is 

A-{x,x) = h{x) + I dtK{t,x), (C.12) 



where h{x) is an arbitrary function of x. The integral can easily be performed in terms of 
multiple polylogarithms. Antisymmetry of A~ {x, x) under an exchange of x and x requires 
h{x) to vanish identically, because 

A~{x,x) = h{x)+ / dtdtA~{t,x) = h{x)+A~{x,x)-A~{x,x) = h{x)+A~{x,x) . (C.13) 

J X 

We thus obtain a unique solution for A^{x, x). 

Having obtained the analytic expressions for A^{x,x) (up to the integration constants 
in A'^{x,x)), we can easily obtain a first-order differential equation for f{x,x), 

dj{x,x) = -[A^{x,x) + A+{x,x)] . (C.14) 

The solution reads 

/(x,x) = j + (C.15) 

The integral can again easily be performed in terms of multiple polylogarithms and the 
antisymmetry oi f{x, x) under an exchange of a; and x again excludes any arbitrary function 
of X only. The solution to eq. flC.14p is however not yet unique, because of the integration 
constants in A'^{x,x), and we are left with three free coefficients of the form, 

(ClC7 + C2C5C2 + C3C4C3)log-. (C.16) 

X 

The free coefficients can be fixed using the requirement that f{x,x) be single- valued (see 
the discussion in Section [7]). Alternatively, they can be fixed by requiring that f{x,x) be 
odd under inversion of (x, x) and vanish aX x = x. We checked that the resulting function 
agrees analytically with the result derived in Section [71 
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